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This is a final report for the AE Guest Investigator Program, NASA 
Grant No. NSG5284, which started July 1, 1978 and ended December 31, 
1982. The research has been under the guidance of Dr. Stanley H. Gross of 
the Electrical Engineering and Computer Science Department of the Poly- 
technic Institute of New York as Principal Investigator and Dr. Carl A. 
Reber of the Goddard Space Flight Center as Co-investigator. 

During the course of the program. Atmosphere Explorer (AE) data 
(AE-C and AE-E) were studied to determine whether coherent wave 
structures were present and to establish their characteristics. These 
data included density, temperature and wind measurements of neutrals, 
and density and temperature measurements of plasma. Power and cross 
spectra have been calculated from the fluctuation data for these various 
geophysical parameters. Both UA and GU data files have been employed. 
Computations were made both at the Goddard Space Flight Center and the 
Polytechnic Institute of New York. It was established that coherent wave- 
like fluctuations are present in all the neutral and ionization parameters 
observed by the satellites, that broad spectra are just about always pre- 
sent extending from 10's km to 1000's km in wavelength in about six years 
of observations. The broad spectra follow power law types of distributions 
with spectral indices different for each parameter. They show that individual 
neutral species are coupled and that plasma is coupled to neutrals. There 
is also a suggestion of flattening of the spectra at the longest wavelengths. 

A statistical study was started to relate these properties to geophysical 
aspects of importance, such as magnetic activity, lime of day, season, 
latitude band, etc. This study was not completed. 

Efforts to find correlative ground based measurements for the 
flight paths of the orbits have been limited and mostly unsuccessful. 

Some ground measurements have been obtained from France of total 
electron content and from the Chatanika radar. These have been partially 
studied and remain to be completed. 

Theoretical models with localized line sources have been developed. 
Numerical calculations have been made for various species in two fluid 
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models. Theoretical cross spectra have been computed for comparison 
with the spectra of measured data. More work is required to improve 
the calculation of cross spectra of the measured data and to model the 
data to obtain possible wave characteristics. However, successful com- 
putation of cross correlation functions have enable the determination of 
transfer functions between geophysical parameter responses. The theo- 
retical models provide insight into the nature of wave properties, re- 
lationships to sources and guidance for interpreting numerical studies 
and measured data. They also serve as test problems for elaborate 
computer programs that are used in models too complex for theoretical 
analysis . 

It was also found that there may be interesting fluctuations of 
neutral parameters where there are "bubbles" in the measured equatorial 
electron density. Furthermore, localized heating at low latitudes over 
restricted longitude bands was found for a number of cases during magnetic 
storms. These findings need further investigation. 

Some of the data imply the presence of quasi- stationary wave 
structures on a global scale. This is concluded from attempts to analyze 
the data with conventional gravity wave theory and from comparison of 
AE-E and AE-C fluctuations and spectra. 

Further work is desirable in terms of (1) completion of the statis- 
tical study of AE-C and AE-E data for levels of fluctuations, spectral 
properties in relation to latitude ranges, time of day, season, magnetic 
activity and other geophysical aspects; (2) the programming of cross 
spectral computations and modeling to interpret the results and to deter- 
mine wave characteristics; (3) the study of French total electron content 
measurements in relation to data from AE-C overflights; (4) the further 
search for corroborative ground based data for comparison with AE 
measurements; (5) the extension of theoretical modeling and numerical 
modeling to guide the interpretation of the data and to relate data to 
source characteristics; (6) the study of the special low latitude heating 
events to relate to the mechanisms whereby energy is transported in the 
thermosphere; (7) the study of "bubbles" and their possible relationship 
to fluctuations in densities, temperature and winds of neutral species; 
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(8) the further study of waves to determine whether a major component 
is quasi* stationary, as implied by some of the data. 


1. The Program 

The program was arranged mainly in three phases: These were: 

1. The analysis of Atmosphere Explorer (AE) data (AE-C 
and AE-E) to determine whether coherent wave structures 
are present and their characteristics. 

2. Theoretical studies of waves emitted by sources in media 
of two or more species with particular emphasis on the 
spectrai region of gravity waves for comparison with 
measured data. 

3. A search for corroborative measurements by ground 
facilities near the tracks of the AE satellites. 


The ground facilities provide temporal data, whereas the satellite pro- 
vides spatial data. Direction of. arrival information may also be provided 
by the former. 

The first two tasks were carried out essentially in parallel. The 
third task was carried out as opportunity permitted. These tasks were 
carried out in the main, but, as typical in research, more work is sug* 
gested to accomplish more fully the original goals. The work, however, 
led to the presentation of material at various meetings of the American 
Geophysical Union, URSI, and at Atmosphere Explorer team meetings. 
Results were presented by Dr. Gross at the NATO/AGARD meeting in 
Naples, Italy, in October 1980. These results were published in the 
proceedings of that meeting. A copy of the paper is attached. Material 
has been prepared for publication in the Journal of Geophysical Research, 
a copy of which is attached as Report No. POLY-MPI- 1427- 83 entitled, 
"Large Scale Waves in the Thermosphere Observed by the AE-C Satellite," 
and another report, POLY- MRI- 1426- 83, has been written covering a 
major portion of the data analysis. This latter report is also attached and 
is entitled, "Gravity Waves in the Thermosphere Observed by the AE 
Satellites. " A number of papers will he forthcoming from this report. 

A number of other efforts will result in additional material to be submitted 
for publication when manuscripts are completed. 
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A great deal of effort was expended on the second task. Both 
analytical and numerical techniques were used. They represent new 
contribution:- to the field, consisting of two fluid models with energy 
sources. Two approaches were used, a theoretical analysis and a 
numerical analysis. The theoretical analysis was to provide insight 
into the nature of the wave solutions and to provide a simpler problem 
to check the computer methods which are intended for more complex 
problems . Further work is required for both. The theoretical analy- 
sis requires a statement of final equations in terms of known functions 
that can be numerically computed. The main features of these formulas 
can be computed and studied to gain physical insight. This work will 
be submitted for publication once completed. 

The numerical analysis of the first version that has been treated 
is nearly complete. It differs from the theoretical analysis by incor- 
porating physical temperature profiles, mean mass variation with altitude 
and more exact formulation, features too difficult to treat analytically. 
Nevertheless, some further calculations are needed for other models 
and for checking with known solutions from the analytical work. It is 
hoped that this work may be completed on receiving additional support. 

The third task has been mainly one of frustration. It was found 
that incoherent backscatter facilities were not operating during particu- 
lar times of good satellite observations. These stations included Arecibo, 
Jicamarca, Millstone Hill, and St. Santin. Ionosonde stations provide 
data that are more difficult to use in a comparison with satellite data, 
though a trio of stations were recently setup by Dr. Millet Morgan of 
Dartmouth College to detect gravity waves. Ui fortunately, these were 
not functioning during most of the period for the AE data. It is possible 
that later AE data (1980, for example) may be used with the Dartmouth 
data. This is an area for consideration in a new program. Ground mea- 
surement data were found as a result of contacts with Dr. Francois 
Bertin of CNET in Paris, France. Measurements were made in France 
of total electron content along a line to a synchronous satellite. These 
were made for several months on the ground in Brittany in 1976 by 
Dr. H. Sizun of CNET. There were overflights of AE-C during this 
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period. The French data were obtained early in 198 1 and analyzed, but 
the comparison remains to be accomplished. The unavailability of the 
Sigma 9 has held up some of this analysis. It can be dc le on another 
computer if a new program is funded. Contacts were also made early 
in 1981 with Drs. R, Vondrak and O. de la Beaujardiere of SRI with re- 
spect to the Chatanika radar data. Some Chatanika data were provided 
at times of overpasses. These data are near the auroral source region, 
and it is important to identify wave disturbances, their nature and charac- 
teristics in this region. Time did not permit further effort here as well, 
and work remains for the future. Other efforts might be fruitful; for 
example, Kwajalein radar data may be useful for comparisons with 
AE-E data. A search for coordinated ground and space measurements 
may bear other fruit. 

There were other discoveries as a result of this work. For ex- 
ample, there is a suggestion that at times of so-called "bubbles" in the 
equatorial region, neutral densities and temperatures fluctuate. This 
requires more investigation. Further, the popular belief that gravity 
waves seed "bubble" phenomena needs clearer demonstration . Studies 
of AE data may provide such evidence. This area is again one for future 
work. Localized heating events at restricted longitude bands at low 
latitudes during magnetically active times were also noted in AE-E data. 
This discovery led to the realization that comparisons for the same 
time of AE-C and AE-E data and with ground based measurements, 
measurements by other satellites may lead to clearer demonstrations 
of the nature of energy transport from high to low latitudes, much more 
than the observation of TID's. This too is an area for continued work. 
Examples of the "bubble connection" and the heating event are given in 
Section 2. The work done on fluctuations, spectra and cross correlations, 
as well as comparisons of AE-E and AE-C data, suggests that much of 
the large scale wave structure seen by the satellites are what can be 
called "quasi- stationary"; that is, waves excited by auroral sources 
which propagate at about 250 km reflect on themselves at the anti- 
podal point. In essence, they are natural spherical modes. There are 
implications here as to loss mechanisms. Further, it is likely that they 
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are generated at F region heights, though this is not certain, from 
particles associated with diffuse aurora, and so may relate to magneto- 
spheric processes, though, again, this is not certain. Section 2 is a 
summary of papers presented at meetings of scientific societies and 
papers prepared for publication. 


2 . Bubble Phenomena and Localized Heating 

Figure 1 is an example of AE-E data for April 28, 1977. It shows 
the variations of densities in the upper graph against time in seconds (UT) 
of oxygen, nitrogen, helium, argon and the electrons. The middle graph 
shows the variations of the neutral temperature and a small portion of 
the ion temperature. The temperature of the MSIS model (Hedin et al. 
1977) is also shown. The lowest graph is for neutral winds of nitrogen 
in the cross-track direction which is mostly in the north- south direction. 
The data are based on 15 second samples and the time span shown is for 
one orbit. Other scales along the abscissa are for magnetic latitude 
(MLAT), dip angle, local solar time (LST), longitude (LONG), latitude 
(LAT), and altitude (ALT). This day was a relatively quiet day mag- 
netically. The mean A index and sum of k indices are shown along with 

P P 

the 10.7 flux. Six vertical, dashed lines are shown in the figure. Four 
of these, those to the left, are at the times of sudden drops in the electron 
density, at the so-called "bubble" features. One sees at these times 
suggestions of corresponding changes in the neutral temperature and the 
wind. Corresponding fluctuations in the neutral densities are too small 
in the figure to be useful. 


One may gain a better sense of the fluctuations of the neutral den- 
sities from Figure 2. These are graphs of the fluctuations of all the 
data of Figure 3, except for the ion temperature which are too sparse. 

All are computed by the filtering techniques described in the attached 
papers. All, except the wind, are normalized to their background 
trends. The bubbles may be observed in the electron density fluctuations 
where there are large negative excursions. Five of eight vertical, dashed 
lines shown in the figure are at the bubbles. Counting from the left, these 
are the third to the seventh line. All except one correspond to the lines 
in Figure 1. The extra line was missed there, but became clearer from 
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the fluctuations. There seem to be perturbations in the neutral densities, 
as well as the temperature and wind, at these times, not necessarily in 
phase, but associated. 

The evidence in Figures 1 and 2 are not conclusive. Much further 
study is needed, particularly of higher resolution data. 

An example of an unusual heating event at equatorial latitudes in a 

restricted longitude range is shown in Figure 3. The day is May 3, 1977. 

It is quite active magnetically. The three hour k indices are shown. 

P 

The previous day was even more active and corresponding graphs appeared 
like those in Figure 3. This figure contains information on density and 
temperature variations as in Figure 1, but there are no wind data and 
there are much more ion temperature data. The event is seen to occur 
centered at the dashed, vertical line on the left in the figure. Note the 
decrease in densities of O, He ana me electrons, and the increase in Ng 
and A. Note the significant increases in both the ion temperature and 
the neutral (N2) temperature. The spread of the event is from about a 
longitude of 136° west to about 60° west. The corresponding magnetic 
latitude spread is from about 16° through 28° to about 25° north. These 
changes are only explainable as due to localized heating, as borne out by 
the temperature variations. Ordinarily, at equatorial latitudes during 
magnetic storms, no particular feature is observed as a function of 
longitude. Only a general warming is apparent (Mayr and Hedin, 1977). 
However, Trinks et al. (1975) from ESR04 data and Trinks et al. (1976) 
from ESR04 and AEROS-A data have suggested similar events though 
these satellites were In highly inclined, elliptical orbits. Other such 
events have been found in the AE-E data. An investigation was starced 
of AE-C data for these days, data from other satellites and some 
ground stations. It is not clear whether the heating is due to a low 
latitude heating phenomena us has been suggested by some, or by a 
"tongue" of precipitation, waves or winds in a restricted longitude 
range extending well down in latitude from the northern auroral 
region. These matters deserve further study. 
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abstract 

We report here on efforts ye 11 lx 1* Atmospheric Explorer (AC) satelllce dace co eacabllsh whether coherent 
waves in Che gravity wave range are praeenc In both neutral and Ionized madia In che thermosphere. AC-C 
dace in particular era shown (altitude 2 50k®, circular orbit inclined 68*). Data consists of In situ argon, 
helium, nitrogen and oxygen densities, plasma density and Ion end electron teoperaturee. Filtering 
provides the fluctuation signals for each which are spectra* analysed for power and cross spectra. The 
observed frequencies ere essentially proportional to che spatial wavenumbers along the satellite track. 
Scale sizes range fro* thousands co tens of kilosacers. 

Modeling efforts to determine the temporal characteristics, direction of propagation and other wave 
parameters from the spectra were made with little success. Unfortunately, efforts to find nesr simul- 
taneous ground measurements co corroborate wave chart ;terlstlcs have been fruitless so far. On reviewing 
the data it was found that turbulence-llke structure over the full range of wavelengths la present with 
different spectral indices for each measured density rad temperature. These appear to be consistent in 
ell date studied to date. Wave packet e ere preeent superimposed on cop of this turbulence structure in 
all data channels, (densities and temperatures) , aany of which exhibit excellent correlation with each 
othar. 


1. INTRODUCTION 

Slnca tha hypothesis of Martyn (1950) 30 yaara ago, it has baan generally accepted chat aany. If not all 
traveling Ionospheric disturbances (TID's) originate from gravity waves in tha neutral atmosphere, lone 
are presumed to be dragged by the neutrals, and In this way exhibit tha wave character! Jtlcs of thi neutrals, 
but with concralnce due to tha earth's magnetic field. Zn spite of aany measurements of meteor trails, 
rocket experiments and radar neutral scattering experiments, chare has bean almost no direct evidence to 
confirm this postulate from ground measurements, because of tha difficulty ot remotely measuring che 
neutral ladium in the upper atmosphere. Tha wide acceptance resulted from Indirect evidence, uch as the 
measured velocity, •cals lengths, periods and the direction of phase propagation (Hines, 1960). 

Direct evidence is more likely in space, elnce spacecraft provide a aultablt platform for making observa- 
tions of both neutrals and ionization. Yet the simultaneous observations of both media has only become 
available In more recent years so that published results are few. Two examples ere Dyson ec al (1970) 
who use Explorer 32 data and Rebar at al (1975) using Atmosphere Explorer C (AE-C) data. In tha first 
reference tha fluctuation in the total neutral density la compared with the fluctuation in the electron 
density. In che second, two examples are given; one compares the fluctuations in the densities of argon, 
nitrogen and helium, and che ocher compares che fluctuations in the densities of ergon, nitrogen and the 
Iona and the fluctuation in tha electron temperature, loth examples were for elliptic orbits for which 
the altitude changes rapidly. These limited cases depend on visual examination of the records to convince 
oneself chat consisted motion exist i. 

Nevertheless, it la apparent In tha examples of the latter paper Chet the fluctuations of helium are quits 
different from those of nierogen or argji in that helium appears at times co be about 180* out of phase 
with the ocher two geese. Also, che amplitudes are different with the amplitude of argon about twice that 
of *2 **d chat of helium about half chat of N2. This behavior has been found co be consistent with che 
theory of gravity waves Involving more chan one species for thermospheric conditions (Grose and Eun, 1976 
and 1978a; Dudis and Rebar, 1977; Dal Canto at al, 1979). 

A more recant exempli involving just neutrals (Hoagy at el, 1979) utilized AE-E date. The interpretation 
and analysis was more sophiaticaced than a visual treatment. In that work a wave is identified and ell 
wave characteristics are derived. However, the theory of their modeling la not fully explained. Further 
examples for neutrals vara given by tha authors of this paper at various meetings of the American 

Geophysical Union and URSI (Groan and Rebar, 1979; Gross at ai, 1979; Huang ec al, 1980) based on A£-C 

data. Soma of chase results will be presented in chit paper along with never material. 

We have been studying data from many orbits of che A£-C satellites, those that are near circular and at an 

altitude of about 250 ka. Circular orbits vers chosen co eliminate the effects of altitude variations on 

the observed f luctuetions. There is evidence of correlation in che fluctuations of all neutral specias 
and ionisation In every case examined. Since the satellite moves very rapidly (• 7.8 km/sae) through ehe 
medium, ehe measurement* essentially provide tha spatial distribution of tha measured quantities. The 
variations exhibit scale lengths typical of gravity waves and art presumed to be such waves. Boeh large 
scale (thousand* of kilometers) and smaller scale (hundreds of kilometers) waves are observed. The scales 
of the waves art determined by spectral analysis using power spectre and cross spectra derived by the 


entropy method (Ulryck and Bishop, 1975). Generally, che stronger spectra are for the larger 
scale eaves. 

The AS satellites provide In situ measurements of the density and drift velocity of argon, helium, nitrogen 
and oxygen as aeasured by the Neutral Atmospheric Composition Experiment (MACS) ; neutral temperature for 
these species from the Neutral Atmospheric Temperature Experiment (NATE); plasma density and ion temper- 
ature by means of the Retarding Potential Analyser (SPA) experiment; plasma density and electron temper- 
ature by means of the Cylindrical Electron Probe, as well as othev experiments (Dalgarno et al, 1973). 

Many of these measurements are essentially simultaneous, making direct comparison of the aeasured 
properties and their fluctuations meaningful. Data utlllsad art available in two different sample 
intervals: 13 sacond averaged data from tha so-called "tJA" files and approximately 1,3 second samples 

from tha so-called "GO” files. 

Our goal was to determloa whether observed wave structures are auroral in origin, so we emphasized A£-C 
circular orbit data becausa that satallite la in a high inclination (- 68*) orbit. Though A£-E has yield- 
ed much data, its low inclination orbit is more suitable for non-auroral-connected events, which are of 
Interest as veil. 

2. EXAMPLES OF DATA AND SPECTRA 

Examples shown in this paper are illustrative. They are taken from data measured during two orbits of 
AJE-C, orbit 3133 on January 20, 1973, and orbit 3339 on February 3, 1973. In both casae tha altitude was 
approximately 250 km. Data for complete orbits are shown baaed upon 13 second averaged data samples, as 
veil as che higher data rate 1.5 second samples for tha study of regions in each of tha hemispheres. 

Complete orbit data are shown in Figure 1. This figure is a composite aede up of the time variation in 
orbit of a number of data channels for orbit 5133, consisting of nitrogen particle density, elsctron 
density and ion and altctron temperatures. Universal time ir seconds is along tha abscissa and ranges 
from 28000 to 34000. The ordlnatee for all the channels are logarithmic* Fluctuations are evident in all 
cases, but correlations between these fluctuations is not clear from this record. (Vertical lines in che 
N 2 and Ng (CEP) graphs are due to spurious computer effects.) 

One obtains the fluctuations by extracting the low frequency or background content from the data shown in 
Figure 1. This Is done by passing the data through a lowpass filter and subtracting its output from the 
original data. Such filtering Is readily done digitally. We utilize the relative fluctuation rather chan 
the fluctuation. The relative fluctuation is obtained by dividing the fluctuation by the background, the 
output of che lowpaas filter as a function of the time in orbit. Figure 2 shows the relative density 
fluctuations for the complete orbit 5359 for nitrogen, argon, helium and the Ionization. The figure also 
provides the fluctuation in che electron and ion temperatures. In contrast vith the variation in total 
quantities in Figure 1, che fluctuations in Figure 2 exhibit, on visual examination, some degree of 
correlation. 

Stronger evidence of correlation may be observed in the power spectra for these same data channels, not 
shown here. These are plots against tha Dopplar frequency In Ez along the abscissa. Since the satellite 
moves so much faster than the phase velocity of a wave, this Doppler spectrum Is essentially kvcosd, 
where k la che wave number, v la che satellite velocity and 9 Is the angle between the wave normal and the 
satellite path. Therefore, for fixed v it rep rr. senes tha component of the spatial dirtributlon as given 
by kcos8. Such spectra are calculated at evenly spaced frequencies by the computer and shown as lines. 
The.c are many common lines in che spectra of che various data channels, particularly for the stronger 
lines. Strong lines are found for the Doppler range from about 0.003 Hz to about 0.008 Hz. For a 
satellite velocity of 7.8 ka/sec this frequency range corresponds to a horizontal wavelength range of 
about 2600 km to 97$ ka. , assuming the satellite track and vavenormal are coincident. Since 9 may not be 
zero, che actual wavelength range may be smaller. Apparently such long wavelength waves are present on a 
worldwide (full orbit) scale. There are also common spectra In the hundreds of kilometer range, chough of 
ouch smaller intensity. 

The shorter wavelengths (hundreds of ka) are much more apparent in the spectra of 1.3 second data samples. 
An example of power spectra for orbit 5133 for the northern latitude range 55* to 0* is shown In 
Figure 3. The Nyqulsc frequency is 1/3 Hz, and the strongest lines are evident from about .025Hz to about 
.040 Hz, corresponding to horizontal wavelengths of about 300 ka down to about 200 ka. Smaller scale 
structure are also evident. 

Further evidence for correlation of the various data channels for orbit 3133 in che saae latitude region 
(35* to 0*) is shown in Figure 4. Here the co- and quadrature spectra are shown for each data channel 
relative co nitrogen. (The power spectrum of a data channel x(t) of total staple time T is 
(^(f) - 1/T| X(f) 1 2 , where X(f) is che Fourier transform and f is the frequency. The cross spejtrum 
between two data channels x(t) and y(t) for total data samples of time T is given by C X y(f)- X (f)Y(f)/T, 
where 'X(f) and Y(f) are the Fourier transforms of x and y and the asterisk means complex conjugate. The 
real part of is che co-spectra, and the imaginary part is the quadrature spectra (Jenkins and Watts, 

1961 ). V It is evident from the strong lines of these spectra that the various channels are not in phase, 
phaseshlft being given by the arctangent of che ratio of the quadrature speccrua to the co-spec trua for 
the satre Doppler frequency. 

The amplitude and phase spectra are derivable from the co- and quadrature spectra, the first by taking* che 
square root of the sua of the squares of the two spectra and the second as explained above to obtain the 
phaseshlft. These spectra are of use in modeling the wave perturbations and their characteristics. 3y 
dividing the amplitude *o*ctrua by the power spectrum of nitrogen at each of che major cross-spectral 
frequencies, che amp I of che relative response of the particular geophysical parameter analyzed 
(density of neutrals or plasma, or esaperature) to that of nierogen is obtained; that Is , | Y(f) /X(f ) \ ,if x 
is che nitrogsu fluctuation and y is the particular parameter of interest. Of course, chi* division only 
makes sense if we are dealing with a coherent wave. The phase spectrua of a channel provides the phase 
response relative to that of nitrogen. For che most significant frequencies, one may deduce che missing 
cnar*:ice?lstlcs of the wave by appropriate modeling using the relative phase and amplitude responses of 



the various channala. Tha Doppler frequency only providaa cha spatial wavelength along tha orbital crack. 

It la not nacaaaarlly cha horizontal wavalangth of tha vava. Thua, by modeling one, in principle , la abla 
Co daduca cha vava 'a temporal frequency, horizontal wavelength, direction of cravel relative to tha 
orbital path, cha velocity of propagation and cha vercleal wavelength. A knowledge of theaa parameters 
together vlch the wave amplitudes la of great laporcance, ainca 1c would permit some daducclona of cha 
energy and momentum c ran* port ad by the vavaa. 

Banda of atrong frequencies are common. Modeling of auch adjacent llnaa should provide redundant informa- 
tion aa though all are from a common source. However, it la a characteristic of the ma xim um entropy 
method (UXry and Bishop, 1975) to cause spreading of these bands and to introduce uncertainty aa to 
their struct Aa a consequence one should not take each line aa exact. Instead, tha average over 

several llnaa may be preferable, (tore work in spectral data analysis in this regard la needad. 

Data ware analysed using cha aodaling techniques as described below. The techniques wars only applied to 
neutrals, because tha thaory of modeling the ionization requires further study incorporating the 
magnetic field* 

3. FIXED LAYER MODELING ANALYSIS 

A full scheme of modeling wave solutions should Include localized sources, multiple constituents and at 
laaat an altitude varying background atmosphert. In addition, winds and horizontal non-uniformity may 
play important roles. A problem of this complexity has not as yet been solved. However, parts of it 
have been treat ad. In the lltaracura solutions have been given for two constituent atmospheres with 
vertical non-uniformity in density (Del Cenio at al, 1979), and for vertical non-uniform density and 
tamperature and mean mass and height dependent collision frequencies (Gross and Eun, 1976, 1978a and b). 
Solutions have also been given for two constituents with all thase non-uniformities and localized sources 
(Gross, 1980). None of these solutions exhibit the degree of variation praaent in the measured data. 

Gross (19 ; ) has shown that the spectra darivad in tha thaoratical probltm for two specias is quite 
smooth anv slowly varying with wavelength and ttmporal frequency. It is not clear then whether the degree 
of variation axt iblcad in spacers aa in Figure 4 is due to the stochastic natura of the source or to 
propagation effects caused by variations in tha madlua. 

Modeling was attempted using a locally hoooganeoua medium about the satellite's altitude to avoid these 
difficulties. We call this oathod tha flxad layer analysis. In this analysis one assumes a gravity 
wave moving through an atmosphere with given background characteristics at an altitude of 250 km: 
temperature, temperature vertical gradient, maan mass, dansity for each constituent and collision fre- 
quencles between species. These values are obtained either from the measured data or alsa from models of 
tha atmos phare for tha particular time of year, location, solar cycle year, etc. The MSIS models (Hadln 
at al, 1977a and b) were used for this purposa. Densitiaa wera obtained by averaging over the measured 
data for the portion of orbit being modeled. Mean mass was also obtained from the average densities or 
alas from the MSX5 models. An atmosphere of more than two species Is treated with two fluid theory by 
taking ona (atomic oxygen), or a combination (atomic oxygan and molecular nitrogen) , as a major species 
and each of tha othars (halium, argon), in turn, as minor. The gravity wave's charactaristics are deter- 
mlned by the major constituent, and the response of helium or argon is ralatable to that of the main 
species from the two fluid theory (Gross and Eun, 1976; Dudls and Rebar, 1977; or Gross, 1980, for example). 
When the main constituent is a mixture of atomic oxygen and nitrogen, tha same two fluid theory provides 
the eigenvalues appropriate to the gravity wave-like solutions for this mixture as well. 

Thus, for orbit 5359 at southern aid-latitudes one obtains the complex relative responses of argon to 
nitrogen and of helium to nitrogen, as explained in the discussion on the cross spectra, for each Doppler 
frequency. These relative responses consist of amplitude and phase information as shown in Figures 5 
and 6, reapectlvely. Here the Doppler frequency is plotted as a harmonic number rather than freauency, 
since cha latter is discrace. Tha fundamental is 1.3021 mHz. The lines in these figures are for the most 
significant responses, all others being weaker, noisy or uncertain due to possible noise. 

With the«e relative responses for each Doppler Frequency, ona has sufficient Information for a model, 
neglecting wind effects and assuming horizontal uniformity, to find the corresponding veve period, 
horizontal wavelength, direction relative to the orbital path and the vertical wavelength. The method of 
defining tha wave consists of scanning an appropriate range of wave periods and horizontal wavelengths to 
find that vava which cornea closest to matching the argon and helium relative responses. A minimum mean 
square error based on the in phase and quadrature components of the responses was used as the criteria; 
that Is, the differences between the calculated and the measured In phase and quadrature components were 
squared and summed for both argon and halium, and Cha sum was minimized. There are two solutions for each 
frequency, ona for an upgoing wave, the second for a downgoing wave. Table I exhibits the wave solutions 
for tha data of Figures 5 and 6 for phaaa-wlae downgoing javes. Wave solutions for 15 Doppler frequencies, 
shown as harmonic numbers , art given. These frequencies are adjacent to each other, and one would expect 
waves of similar characteristics. Instead, tha table shove a comp lata scatter in the periods and 
horizontal wavelengths, (shown as > x ) . Furthermore, one cannot determine the wave direction relative to 
the satellite path, 9, for moat cases, since the best solutions give cos 9 > 1. The measured and calcu- 
lated amplitudes for the argon /nitrogen and hellum/nltrogtn responses, shown as A/N 2 and H«/N 2 in the 
table, art considerably ilfferent. Tha fit is poor; the errors are comparable to the measured values. 
Similar results are obtained for phase-wise upgoing waves, for Che ocher parts of the same orbit and for 
ocher orbits. If on# averages adjacent lines to compensate for the peculiarities of Che maximum entropy 
method and models the resulting values, similar results are obtained averaging three at a time, four at a 
time, etc. Varying the temperature or the collision frequencies by + 10X yields no improvement. 

Thus, the modeling attempts were unsuccessful. The medium may be far too complex for the assumptions. 

Data for chouaands of kllooaters must be uaed to analyze waves of several hundred kilometers wavelength. 
Undoubtedly, winds, wind shears and horizontal non-uniformity are not inconsequential over such distances. 
Further, the greater the distance from the source, the more is the change In wave characteristics 
(Chinonaa and Hines, 1970). To model all these effects requires information not readily available, and, 
even if available, difficult to incorporate. We cannot rule out noise in the data, aa well. 
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When we started thla program ve vara optimistic chac va could determine cha characteristics of the vavaa 
from cha satellite measurements because of cha nuabar of simultaneously aaaaurad quantities. Thasa negative 
results, however, lead ua co ballava chat 1c la essential co corralaca spacecraft and ground aaaauraaanca 
Co find all cha charactarladca of cha vavaa, caaporal and spatial. Unf ortunataiy , co daea va hava aoc 
baan abla co find aultabla ground axparlaancs functioning during a aatalllca paaa ovar cha scadon vhan 
avanca occurred. 

Another approach haa baan co ravlav cha data In a dlffarant faahlon. thla haa lad co cha discovery of 
nav t possibly important, characcarlaclca Chat had not originally baan appreciated, nor poaalbly known ac 
all before, aa far aa ehaaa auchora can da t amine. Though chaaa new raaulca ara preliminary, ehay may 
hopefully laad co a ouch baccar understanding of cha mechanisms at work and chair relationships co 
important geophysical aspects, such a a magnetic aeonaa, solar cycla variations, saaaon, etc. These find- 

in,* will now b« •splninad. ORIGINAL PAGE ES 

4. USGE SCALE HAVE DISTRIBUTIONS AND TURBULENCE OF POOR QUALITY 

By a vary simple process, chac la, by raploecing cha spectral daca shown In Flguras 3 or 4 on a log- log 
graph, an unappreciated faatura became evident. In this paper va show this only for eha power spacers. 

All are for cha northern latitude range 53* co 0* of orble 5133. Va see similar faacuras in othar pares 
of the same orble and for ocher orbits. 

The power spacers of Sj, 0 and alectron density n ft are plotted versus the harmonic frequency number in 
Figure 7 through 9. They are shown aa docs on the graphs. Thoaa below about harmonic numbs r 20 era la 
the lovpsss filter range and should be Ignored; Che fluctuations of lncarasc are chose froa harmonic 
number 20 and up. Harmonic number 20 corresponds co a wavelength of 300 km along the orbital path. Also 
shown on each graph is a base logarithmic acralghc line fit. Such a fit la expressible as An, where A 
Is an amplitude constant, n la the harmonic number and a la the spectral Index characterizing the wave 
number decay race. The fits for Figures 7-9 ere based on the spectrum from n * 25 to n * 100. The 
deviation of tha spectra In these figures la, except for some discrete values of a, relatively smell. 

The values of a for theae three fits are -5.3719 for nlerogen, -3.14128 for atomic oxygen and -2.1401 for 
the electron density. Curiously, the ratio of the value of a for nitrogen to chat for oxygen Is very close 
to the ratio of their molecular weights. Not enough daca have been analyzed to establish whether this is 
real or a coincidence. 

Do the fits in Flgurss 7-9 Imply a turbulence- like spectrum for these large scale waves? The graphs 
arc really wave number plots, and aa such they are reminiscent of eddy distributions (Monin and Yagloa, 
1971), but for much larger scale lengths. It la also possible that such spectra are representative of a 
stochastic sourca and thertfore not turbulence. It may originate froa non-linear Interactions In the 
region In and about the source. In any event, it Is believed chac strong non-linear Interactions muse be 
involved, whether turbulence or source derived. 

Plots similar to those of Figures 7 through 9 were made for helium and cha electron temperature. These are 
not shown here. The fit for helium for the spectrum between n ■ 25 and n * 100, a a in Figures 7-9, haa 
the exponent a - -0.1645. Thla line yields a power spectral density almost independent of wavenumber, 
like white nolee, and corroborates a noise- like distribution evident in data, particularly beyond n • 50 . 

If a fit la made between a ■ 25 and n ■ 50, one obtains a second line with an exponent a - -1.45. A close 
examination of the data shows that there la a decreasing spectrum in this range of a. Apparently at 
shorter wavelengths noise dominates over the decreasing spectrum suggested by the second fit. This latter 
fit Is probably not very accurate due to the influence of noise. Though not shown here as veil, the 
fluctuation va. time data for helium In this portion of the orbit, unlike Chac in Figure 2 for helium, 
appears quite noisy, implying that the MACE instrument was limited by the noise level in this measurement 
of helium. Nitrogen and oxygen are more abundant, so that noise does not appear co be a problem. The 
electron density measurements are from another instrument, the cylindrical electron probe. 

The spectrum for the electron temperature is even noisier chan that of helium. Two fits are also suggested. 
One Is for n • 25 to n • 50 for which the exponent is positive and small (a • 0.342), as is expected for 
white noise. The other fit is for n - 25 to n • 40. this fit was computed because there seems to be a 
decreasing spectrum in the range. Its exponent Is a * -2.31. What is interesting is that the elimination 
of juat 10 frequencies froa n ■ 41 co n • 50 products this change. In contrast, chs computation of fits 
for ranges other then n » 25 to 100 for nitrogen, oxygen and the electron density causes minor modification 
of their exponents. As for helium, the correct fit for the electron temperature spectrum is uncertain* 

Ac best, one can only conclude there is some hint of a dscreaalng spectrum that becomes overwhelmed by 
noise as the scale size decreases. Aa far as chess authors know, such an electron temperature spectrum is 
the first of its kind ever deduced in natural plasma. Unfortunately, room does not permit Its inclusion. 

Over the range n - 23 to n * 50 thare appear to be significant departures in the various spectra from their 

straight line fit ac or near particular values of n. These departures or deviations may imply a wave 

system sitting on a turbulence-like spectrum. The features are made evident by plotting the difference, 
or residuals, vs. n on linear plots. Flgurs 10 snows these differences for nitrogen, oxygen and the 

electron density. The match of the peaks in these differences, between those of nitrogen and oxygen, is 

excellent, and that for the electron density is quite good. We find this match even though the maximum 
entropy method tends to broaden any peaks. The first large peak in all three of the residuals is at a 
wavelength along th r, orbit of about 272 km. The second peak (missing In the electron density spectrum) 
is at about 240 km. The third, fourth and fifth peaks are ac about 190ka, 150ka f and 133 km, respectively. 
Though speculative, this group could be Indicative of a wave packet moving through the medium. 

The deviations for helium and the electron temperature are shown in Figure 11, based upon the n • 25 to 50 
fit for helium and the n - 25 to 40 fit for the electron temperature. Though there is some degree of 
overlap with the jpectra in Figure 10, the wide excursions are not surprising in view of the noisiness of 
the data for these two channels. These spectra are not considered too reliable. 



3 . 


CONCLUSIONS 


The a«|«clvi results chat were obtained in attempting to aodal tha fluctuations in cha neutral a aa gravity 
wavaa la not surprising in vlav of tha complexity of cha vava structura saan In Flguraa 7 to il. 

Apparantly, non-linaar processes, probably aaaoelatad with strong stochaatlc sourcaa in tha auroral region, 
and possibly from meteorological sourcaa, ara at *ork producing a turbulence- like spectrum in tha gravity 
vava ranga stratchlng from vavalangths of aany hundrada of kilometers down to tha smallest scalaa 
obaarvabla bafora nolsa dominates, on tha ordar of tana of kilometers. Work with coop lata orbit data also 
auggaata that tha turbulence- like structura extends to vavas of thouaand of kilonatars in wavelength. Tha 
axponanta of tha spacers for asch of tha nautral constituents, tha plasaa and alactron temperature ara of 
different valuaa in tha llnitad axanplaa givan hara and tha fav ochars that hava baan invast igatad. Why 
this is cha caaa is not clear, but night ariaa froa non-linaar cascading procaaaaa and loss nachanlsna 
aaaoelatad with diffusion. Thara la much theory to ba invastlgacad hara. It la not claar how or if thasa 
waves ara propagating* Ara tha waves stationary, ara they propagating in all directions or from tha 
sources outwards? Tha deviations of Figures 10 and 11 suggest chat wave packets ara propagating on top 
of turbulence. It la conceivable that such wave packets would be observed as a traveling ionospheric 
disturbance on tha ground. 

A moat interesting future investigation will consist of the statistical study of these characteristics and 
how they may relate to geophysical aspects, such as tha year in tha solar cycle, seasonal effects, 
magnetic storms, solar sector boundary crossings, ate. There is also a need to study similar character- 
istics in the equatorial ragion, as is possible with AE-E data, for relationships to tha equatorial 
electro jet and meteorological events. 
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2.6 

18.8 
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Figure 1 Densities and temperatures 
for one entire orbit, 
orbit 5133. 


Figure 2 Relative fluctuations of 

densities and temperatures 
for orbit 5359. 
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Figure 3 Power spectra for a part of orbit 
5133, 55* to 0* latitude* 


15133 MSM/GU 29600 


CO -SPECTRA 



FREQUENCY 


QUADRATURE 

SPECTRA 



0667 2000 3333 

FREQUENCY 


Figure 4 Go- and quadrature spectra for the 
data channels of Figure 3 relative 
to nitrogen. 
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Figure S Amplitude ratios of argon to nltrogan 
and ha Hum to nitrogen for significant 
lines of the cross spectra, orbit 5359, 
southern mid-latitudes. 




Figure 6 Phase spectra associated 
with the amplitude 
spectra of Figure 5. 
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Figure 7 Power spectrum of nitro- 
gen vs. tne harmonic 
number n for data of 
Figure 3. Best fit line 
based on the spectra 
from n • 25 to 100 is 
shewn. 


Figure 3 The same as 
Figure 7 but 
for atomic 
oxygen . 


Figure 9 The same as 
Figure 7 but 
for the 
electron 
density. 
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Figure LO Plots of residuals of spectra Figure il The same as Figure 10 but for helium 

in Figure 7 - 9 in the range and the electron temperature, 

n ■ 20 to n ■ 50. 
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Atmosphere Explorer (AE) satellite data were used to Investigate the 
spectral characteristics of wave-like structure observed In the neutral and 
Ionized components of the thermosphere. Power spectral analysis derived by 
the maximum entropy method Indicate the existence of a broad spectrum of scale 
sizes for the fluctuations ranging from tens to thousands of kilometers. A 
new discovery is a power law decrease of power with decreasing wavelength, 
much like a turbulence structure, though not necessarily turbulence, which is 
generally observed in all data examined. Waves that may be gravity waves are 
observed in the species and other parameters studied. These appear superposed 
on top of the turbulence-like structure. Correlation of many waves in the 
spectra of the various parameters, including neutral and plasma parameters, is 
quite evident, demonstrating strong coupling between the constituents of the 
medium. A typical example is shown which illustrates these features. 
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The direct observation of large scale waves both In neutral and Ionized 
media In the thermosphere has been made possible by the measurements of the 
Atmosphere Explorer (AE) satellites. These satellites provide in situ 
measurements of the density of argon, helium, nitrogen, and oxygen, the 
temperature and drift velocities of nitrogen or oxygen, plasma density, ion 
and electron temperatures, as well as other measurements (Dalgamo et al., 
1973; Pelz, et al., 1973; Nler et al.- 1973; Spencer, et al., 1973; Brace et 
al., 1973; Hanson et al., 1973). Many of these measurements are essentially 
simultaneous, making direct comparison possible. Data from a number of orbits 
of the AE-C satellite have been studied, iocu 3 ing on those that are near 
circular at an altitude of about 250 km. Circular orbits eliminate the 
effects of altitude variations on observed fluctuations. The AE-C orbit is 
Inclined *68° . 

There is evidence of correlation in the fluctuations of all neutral 
spec! s and ionization in every case examined, implying the existence of wave 
like structures, aud supporting previous, more limited examples based just o a 
visual examination of the measured data (Reber et al (1975) using AE-C data: 
Dyson et al (1970) using Explorer 32), and analytically treated examples 
(Hoegy et al 1979). Since the satellite moves very rapidly (— 7.8 lon/sec) , 
the measurements essentially provide the spatial distribution of the measured 
quantities. The scales of the waves are determined from spectral analysis 
using power spectra derived by the maximum entropy method, MEM (Dlryck and 
Bishop, 1975). Scale lengths are observed from 10's to 1000's of kilometers. 



Data and Spectra 
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One typical example of AE~C data is shown In Figure 1* These data are 
approximately 1.5 second samples, plotted vs time over a time Interval of 750 
seconds. They are from orbit 5133 on January 20, 1975 for a portion of the 
orbit in the northern hemisphere from about 10° to 54° latitude. Altitude 
varies from about 244 las to 252 lan, and the local time ranges from 1100 to 
1300 hours. Magnetic and solar activities at the time were relatively 
quiet. The indices were Ap » 7, Kp • 2, AE ■ 114, and the 10.7 cm solar flux 
Fjq # 7 *77. Shown are the densities of oxygen, nitrogen, helium and the 
electrons and the electron temperature plotted against time in seconds (UT). 
The ordinate scales for the densities are logarithmic. 

Some fluctuations are evident in the geopnysical parameters in Figure 
1. With the exception of nitrogen and oxygen, and possibly the electron 
density, any relationship between these fluctuations Is not clear from the 
figure. Relationships become more evident on separating the fluctuations from 
the data. This is done by extracting their background, or low frequency 
content using a low pass digital filter. Its output is subtracted from the 
data to obtain the fluctuations. Since the scale of the fluctuations varies 
with the parrameters, comparisons between parameters are best made by 
normalizing the fluctuations. This is done by dividing the fluctuations by 
the output of tne low pass filter. The normalized fluctuations (designated 
An/n for densities, or AT/T for temperature) for each of the parameters in 
Figure. 1 are shown in Figure 2. 
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. The filter used is a Hamming non-recursive filter including Lanczos sigma 
factors (Hamming, 1977), Its order is 60. It passes all frequencies below 
0.01 Hz, equivalent to 780 km wavelength. The attenuation of the filter at 
0.01 Rz is 0.5. Less than 1Z of the power Is passed at 0.02 Hz and far less 
for higher frequencies. Finer details of its transfer function are not given 
here. 

The power spectra for each of the fluctuations in Figure 2 are shown in 
Figure 3. These are plots of power spectral density as a function of 
frequency as observed by a moving satellite (analogous to a Doppler 
frequency). About 380 fluctuation samples for each parameter are used after 
filtering to obtain the spectra by means of the MEM with an order equal to 100 
and using the Burg algorithm. This order was selected after study of the 
effects of the MEM order on spectra and by comparing with FFT spectra and 
spectra derived by other means. The details are not given here. The spectra 
are shown at evenly spaced frequencies in the figure, making them have the 
artificial appearance of line spectra. The frequency range shown extends o .t 
to the Nyquist frequency, 1/3 Hz. A wavelength scale is shown below the 
frequency scale. This scale corresponds to the equivalent wavelength of a 
wave standing along the satellite track and is obtained by dividing the 
satellite velocity, 7.8 km/sec, by the corresponding frequency. A wave of 
given wavelength traveling at an angle to the satellite track would appear to 
have a longer wavelength along the track. 

More features of the spectra are brought out by plotting on log-log 



graphs. Figure 4 shows the power spectra of oxygen, nitrogen and electron 
density replotted in this fashion. The crosses are at frequencies that are 
multiples of 1.302 mHz, which corresponds to the time period, 768 seconds, of 
512 samples spaced 1.5 seconds apart. The frequencies shown are mostly above 
the filter range, so that these spectra represent the high frequency content 
of the fluctuations. The corresponding wavelengths of the spectra shown there 
are from about 545 km to about 40 km. Shown on each graph are straight lines 
which represent least squares power law fits. These fits are based on the 
spectra above the filter range (for wavelengths below 374 lan) . The exponents, 
or spectral indices, differ for all three species shown: -5.38 for nitrogen, 

-3.51 for oxygen, and -2.25 for the electron density. 

The spectra of helium and the electron temperature were not included in 
Figure 4 because of their noisy appearance. They exhibit many of the features 
of white noise, particularly at higher frequencies. Nevertheless, their lower 
frequency spectra have some features in common with those of oxygen, nitrogen 
and the electron density, in spite of the noise, as may be seen in figure 3, 
but the details are not given here. 
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One observes In Figure 1 relatively smooth variation in the densities of 
N 2 and 0 and rather irregular, or noisy, variations in the density of He. The 
noisiness of the He data is well brought out by the appearance of its 
fluctuations, in contrast with those of N 2 and 0, as shown in Figure 2. It is 
also seen In the broad, white-Qoise-like spectrum of He, In contrast with the 
narrower band spectra of N 2 and 0, as seen in Figure 3. Nevertheless, the 
lowest frequency peaks in the He spectrum are about the same as those of N 2 
and 0 and one can associate some of the fluctuations of He in Figure 2 with 
those of N 2 and 0. They appear out of phase. 

The variation of the electron density in Figure 1 also exhibits 
fluctuations along with large scale changes, of which the largest, at about 
30250sec (UT), is apparently associated with the equatorial anomaly. There 
are also sections of the data that exhibit fluctuations much like those of N 2 
and 0. The relative fluctuations of the plasma density in Figure 2 appear to 
be a bit noisy, yet somewhat like those of N 2 and 0, as may be seen in that 
figur . Figure 3 shows that the power spectrum of the electron density is 
also limited In width, similar to N 2 and 0. 

Considerable variation in the electron temperature is also apparent in 
Figure 1. Rapid, noise-like fluctuations are evident, as well as a large 
scale change at about 30250 sec (UT) that is associated with the anomaly. Its 
fluctuations, as seen in Figure 2, also appear noisy, though perhaps not as 
much as that of He. One can associate some of the fluctuations with those of 
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N 2 and 0. Its power spetrum In Figure 3 Is also broad, like white noise, 
though somewhat reduced at higher frequencies. The lower frequency peaks are 
about the same as those of Np and 0. 

The fluctuations of 0, Np, and the electron density exhibit wave-like 
structures, as seen in Figure 2. N 2 and 0 fluctuations are quite similar, 
though of different magnitude. One observes structure with a number of 
periods, but there is an outstanding periodicity roughly of the order of 35 
seconds, or about 270 kilometers on converting to the corresponding path 
length along the satellite track. The level of the N£ fluctuations averages 
about 3Z with some fluctuations reaching peaks of as much as 7 to 8Z. The 0 
fluctuations average about 1Z with some peaks equal to about 3Z. These 
fluctuations are not particularly intense, though, the example is typical and 
adequate for illustration here. The electron density fluctuations are larger 
than those of N£ and 0, averaging about SZ with some peaks reaching 13Z. The 
noise-like fluctuations of He and the electron temperature both average on the 
order of 5Z. A close examination of the fluctuations of Np, 0 and the 
electron density shows very small phase shifts between their various wavelet6 
with the electron density leading by about 20 seconds. 

The power spectra in Figure 3 of Np, 0 and the electron density exhibit 
peaks at about the same frequencies, or wavelengths, corroborating the 
similarities in the structure of their fluctuations. The first spectral peaks 
are at frequencies on the order of 0.01 Hz, within the low pass filter band 
and therefore attenuated. Their wavelengths are on the order of 800km. A 
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larger second peak, a still larger third peak and a much smaller fourth peak 
are evident in all three spectra. These are at wavelengths of roughly 430, 

270 and 195 km respectively. The 270 k» wavelength corresponds to the 35 
second periodicity seen in the fluctuation data. It is believed that all 
these peaks are due to gravity waves, though it is not clear if they are 
propagating, and, if so, with what velocity, direction and temporal 
frequency. Such wave characteristics are more easily determined in 
conjunction with ground measurements which, unfortunately, were not available 
during passes of the satellite over such facilities. It tnay be possible to 
obtain this information from cross spectra (Hoegy et al., 1979), but 
difficulties with the MEM algorithm for obtaining such spectra has precluded 
their use for the purposes of this paper (Gross et al., 1980). 

The log-log presentations of the power spectra of N 2 , 0 and the electron 
density in Figure 4 indicates that these fluctuations consist of broad 
spectra. The power law fits illustrate that the power in each specie 
decreases with increasing frequency, or decreasing wavelength. The decrease 
of power with frequency is greatest for nitrogen, less for oxygen and least 
for the electron density. We believe that part of the power law - like 
variation is due to the sources. Why the rates for each are different is 
unknown. If these spectra arise from the same sources, as seems likely, it is 
conceivable that the differences arise from the propagation characteristics of 
the medium. Such power law distributions are also reminiscent of eddies or 
turbulence (Monin and Yaglom, 1971), though of much larger scale lengths, 
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leading to the description "turbulence-like". However, not enough information 
la available to establish whether or not turbulence is Involved. It is 
possible, for example, to obtain similar spectra from simple perturbations. 
Since coherence tends to narrow spectra, the breadth of the spectra may merely 
be indicative of non-coherent sources. 

Peaks are evident in the discrete spectra in Figure 4. Those 
corresponding to wavelengths greater than 100 km are listed in the table 
below. All figures are in kilometers. 



Though the values are to three figures, these peaks approximate the 
wavelengths of the waves along the satellite track. This is due to the use of 
discrete frequencies and because of offsets arising from the MEM. The longest 
wavelength is in the filtered range and therefore somewhat attenuated. All 
others are unattenuated. The first three correspond to peaks clearly seen in 
Pigure 3. The match between these various wavelengths for the three species 
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is excellent, particularly for nitrogn and oxygen* (Though not given, the 
corresponding peaks of He and the electron temperature also match well*) 

Those of the electron density are either the same or close, with the exception 
of a missing (or possibly indistinct) peak at about 130 km. These spectra 

provide evidence that common wave systems are affecting all three species, and 

probably the other parameters as well. Also, this example clearly 
demonstrates the coupled behavior of the ionosphere and neutral atmosphere. 

Though expected, such evidence is generally sparse. It is likely that the 

waves in the ionization result from wave motion of the neutrals. All these 
waves may be medium scale gravity waves, as described by Georges (1968). They 
appear to be superposed on the turbulence-like spectra of each of the species. 

More than one source may be needed to produce the medium scale waves; for 
example, one source for each wave. Each source would excite its own wave 
characterized by its own structural size. Multiple sources may occur 
typically during intense meteorological activity. On the other hand, if the 
waves emanate from one source they might represent different structural 
features of the source. Regardless of the nature of the sour cer, waves of 
these scales are likely to be local, since waves of wavelengths of only 
several hundred kilometers or less would be well attenuated, if they 
propagated from distant sources. 

The spectra shown here are just for one example from AE-C. We find 
similar spectra for other parts of the same orbit and for other orbits in all 
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the data we have analyzed. The length of Che data samples In Figure 1 is 750 
seconds. This is about the time Interval (representing about 6000 kilometers) 
required to study waves for specific latitude belts of geophysical 
significance, such as mid-, high and equatorial. Scale lengths on the order 
of tens to hundreds of kilometers are found. When longer time samples are 
analyzed, the scale sizes are larger. Full orbit data provide similar 
spectra, but sizes to thousands of kilometers are obtained still characterized 
by power lav fits and peaking at Individual frequencies as in Figure 4. Such 
spectra, however, represent a mix of '’local" phenomena as well as worldwide 
large scale phenomena. A broad, power law-type spectrum with waves 
superposed, a new discovery for waves of these scales, is apparently 
omnipresent as far as we can tell from the data analyzed. 

The nature of the spectra described here are also corroborated by another 
experiment. Total electron content was measured between two spacecrafts 
during the Apollo-Soyuz Test Program by the differential doppler technique 
(Estes and Grossi, 1980; Weiffenbach et a!., 1977). These measurements 
exhibit power spectra with similar features with the measured data processed 
directly by the MEM without filtering (Gross et al., 1981). It is also 
Interesting that Dyson et al (1974) found similar structures in AE-C electron 
density measurements. 

The complexity of wave structures is further exhibited by the broad 
frequency spectrum found by Titheridge (1971) from ground temporal 
measurements of total electron content. 
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The material shown here provides further insight into the nature of waves 
in the upper atmosphere as observed by satellites. The data of the A£ 
satellites, both AE-C and AE-E, are rich with such evidence. Much more work 
is needed in analyzing these data to establish the characteristics of these 
waves, particularly gravity waves, their function in transporting and 
distributing energy and momentum, and their properties in terms of geophysical 
parameters, such as solar and magnetic activity, season and time of day. 
Acknowledgement : These authors thank K. Coakley of Computer Sciences 

Corporation for his aid in calculating some of the material in this paper. 
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Figure Captions 


Densities of nitrogen, oxygen, helium and the electrons 
and the electron temperature vs time in seconds, altitude 
and latitude from- January 20, 1975 Atmosphere Explorer 
C Orbit 5133. Neutral density measurements are from 
Neutral Atmosphere Temperature Experiment (NATE), and 
electron density and temperature from Cylindrical Electro- 
static Probe (CEP). 


Relative fluctuations in the densities, An/n, of nitrogen, 
oxygen, helium and electrons, in election temperature, 
AT/T, vs time in seconds for data in Figure 1. 


Discrete power spectra of relative fluctuations in 
Figure 2 vs frequency in mH.. for the entire range to 
333 mHz, the Nyquist frequency, and vs wavelength in 
kilometers . 


Discrete power spectra of Figure 3 from 14.3 mHz to 
195 mHz replotted with logarithmic axes for nitrogen, 
oxygen and electron densities vs frequency and wave- 
length, spaced logarithmically. Power law fits are 
shown in the form f” , where f is the frequency, 
based on spectra between 20.8 mHz and 130.2 mHz. 
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ABSTRACT 


Atmospheric Explorer C <AE-C> data are 
analyzed to study wavelike perturbations in the 
thermosphere at an altitude o-f about 260 km. The 
data are -from a large portion o-f an orbit on Janu- 
ary 20, 1975. The examples shown are typical o-f 
many other orbits o-f satellites AE-C and AE-E. 
Four geophysical parameters are analyzed. These 
are nitrogen and oxygen densities, electron den- 
sity and ion temperature as measured by three dif- 
ferent instruments on board AE-C. 

The data are processed by normalizing them to 
their average values and their fluctuations are 
obtained by passing the normalized data through a 
high pass filter. These fluctuations exhibit 
strong periodicities that persist throughout the 
records, regardless of the filter cutoff fre- 
quency. The fluctuations are cross correlated by 
direct comparison, by comparison of their spectral 
content as obtained by the MEM and fft, and by 
computing cross correlation functions. 

The ratio of the fluctuations and their rela- 
tive phase shifts are treated as the magnitude and 
phase shift of a transfer function between the 
responses of all parameters relative to the 

. response of the oxygen density. These transfer 
functions are used in a plane wave modeling study 
..Ib-deduce all wave characteristics. 1 1 . was not 



possible to account -for the wave characteristics, 
suggesting that the fluctuations may be a mixture 
of standing and traveling waves. ' 

It is found that there are large scale struc- 
tures in all four parameters with horizontal scale 
sizes ranging from about 4000 km to about 400 km. 
Both ionization and neutrals are driven by these 
waves. Ionization, including even the ion tem- 
perature, is coupled to the neutrals for most 
waves. The coupling and wave system persists over 
large global distances and in both hemispheres. 

The spectra for all parameters contain peaks 
at wavelengths that are confirmed by the periods 
of the fluctuations and decrease with decreasing 
wavelength with a power law type of variation. 
The 'spectral index', or exponent of the power 
law, for the four parameters differ, though those 
of oxygen and nitrogen are close in value. Com- 
parison of the spectral peaks show that many waves 
are common to all parameters. 


L. In troduc t i on 

Wave-like per turbat i ons of constituents in the thermo- 
sphere are detected by measurements of the ionosphere by 
ground facilities and by instrumentation on board satellites 
measuring both neutral and ionospheric parameters. Pertur- 
bations of ionization are known as traveling ionospheric 
disturbances <TID's) (Hines, 1974; Georges, 1968; Yeh and 
Liu, 1974). A number of reports of large scale TID's (from 
hundreds to thousands of km) have shown that they propagate 
from high latitude regions during auroral substorms (Thome, 
1968; Davis and da Rosa, 1969; Testud et al,1975>. These 
waves are believed to be caused by heat and momentum input 
into the atmosphere in the auroral zones from particle pre- 
c i p i t i on and Joule heating (Vickrey et al,1982; Wickwar et 
al , 1975). Their origins have been reviewed by Hunsucker 
(1982). Studies have also been made of these waves in asso- 
ciation with magnetic storms by Richmond and Matsushita 
(1975), Mayr andVolland (1973 and 1974), Testud et al 
(1975), Roble et al (1978), Richmond (1979) and Prolss 
(1982) . 

Medium scale TID's (from tens to hundreds of km) are 
also observed, but their sources have not been well identi- 
fied, though it is believed that meteorological phenomena 
may produce most of them. Larsen et al (1982) have associ- 
ated thunderstorms with these waves. Mastrantonio et al 
(1976) have also associated these waves with jet streams in 
the atmosphere, and Bertin et al (1977) have detected waves 
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related to the jet stream. '*• 1 

The characteristics of all these TID's are 'found to be 
generally in accordance with gravity wave theory 
(Hines, 1974) . It is believed that the neutrals are per- 
turbed by gravity waves as a result of localized sources of 
heating and momentum, and the neutrals, in turn, drive the 
ionization by collisions to produce these TID / s. Gravity 
waves are believed to be associated with Spread F in the 
ionosphere and ionospheric bubbles in the equatorial region, 
affecting communications with synchronous sattelites (Ander- 
son et al , 1982; Kelley et al , 1981). Anderson et al*(1982) 
also pose another mechanism for generation of gravity waves 
in the thermosphere in the equatorial zone. 

Gravity wave theory is concerned with two types of 
waves called gravity waves and acoust i c-grav i ty waves. 
Acoust i c-grav i ty waves make up the upper frequency branch 
and gravity waves the lower frequency branch of the hydro- 
dynamical system in which both gravity (buoyancy) and 
compressibility act as restoring forces. The upper branch 
becomes identical to ordinary sound waves in the audio fre- 
quency range, but at its lowest frequencies such waves are 
called / infrasonic waves'. This branch has a lower cutoff 
frequency called the 'acoustic cutoff frequency' , and there 
is a cutoff region below this frequency, which separates the 
two branches. The lower end of the cutoff band, which is 
the upper limit of the gravity wave branch, is called the 
'buoyancy f requency' . The cutoff band is determined by the 
temperature and its vertical gradient and by the mean molec- 
ular mass of the medium, so that it varies somewhat with 
altitude. The cutoff band is relatively narrow, and periods 
within it range from about several minutes at low altitudes 
to from 10 to 15 minutes at higher altitudes, depending on 
temperature. The periods of gravity waves then range from 
the period of the buoyancy frequency to many hours, whereas 
those of acoustic-gravity waves range from the acoustic cut- 
off period to periods of seconds and less. In the gravity 
wave branch, long period waves usually have large horizontal 
wavelengths and shorter periods smaller wavelengths (Hun- 
sucker , 1 982) . 

Though fluctuations in ionization in the thermosphere 
have been measured both on the ground and by satell i te borne 
instrumentation, the detection of wave-like per turbat i ons of 
neutral constituents in the thermosphere depends mostly on 
access by rockets and satellites. Radar is used only up to 
an altitude of about 100 km. Satellite borne experiments 
have detected wave-like per turbat i ons , both in total neutral 
density and temperature; for example, ESRO 4 (Prolss and Von 
Zahn, 1974 a and b; Trinks and Mayr , 1976), OGO 6 (Reber and 
Hedin, 1974; Taeusch et al , 1971), AEROS-A (Chandra and 
Spencer, 1975 and 1976; Trinks et al , 1976), AE-C (Reber et 
al , 1975; Potter et al , 1976). The density variation of 
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individual species, namely nitrogen, oxygen, helium and 
argon, have also been measured on board ESRO 4, AEROS-A and 
AE-C. Vertical velocity variations of these species have 
also been measured by AE-E. 

Until the advent of satellites there has been little 
confirmation that wavelike per turbat i ons in both the neu- 
trals and ionization are correlated. Visual examination of 
data measured on board Explorer 32 (Dyson et al , 1970 and 
1974) and AE-C (Reber et al , 1975) provided some support for 
such correlation. Some additional evidence were provided by 
Hoegy et al (1979) employing some analytical modeling. 
Gross et al (1980 and 1981), found still more evidence on 
the basis of visual comparisons of spectra and fluctuations. 

Here we report on a study of some AE-C data to examine 
these wavelike per turbat i ons in greater detail using various 
data processing techniques. Since instrumentation on board 
AE-C make in situ measurements of many parameters nearly 
simultaneously, it is possible to cross correlate their 
fluctuations to seek their common properties and to show the 
degree to which they are truly correlated. AE-C measure- 
ments are made of neutral densities, electron density and 
temperature, and ion temperature, as well as many other 
parameters. Here we use nitrogen and oxygen densities as 
det ermined by the Neutral Atmosphere Temperature Instrument 
(Spencer et al , 1973), electron density from the Cylindrical 
Electrostatic Probe experiment (Brace et al , 1973) and ion 
temperature from the Retarding Potential Analyzer experiment 
(Hanson et al , 1973). The satellite was in a near circular 
orbit at an altitude of about 260 km while measuring the 
data treated here. The orbit was inclined 68°, 

The satellited velocity in orbit is about 7.8 km/sec. 
The waves of interest have group velocities of the order o x 
the acoustic velocity or less, about 800 m/sec in the ther- 
mosphere. Since the satellite velocity is so much greater 
than the wave velocity, the in situ measurements are essen- 
tially a snapshot of the distribution of the measured param- 
eters at the altitude of the satellite at a particular time. 
Thus, wavelike variations in the measured time series are 
really spatial variations. One observes the horizontal 
wavelength for a circular orbit and not the wave period, 
though the terns frequency and periodicity will be used 
often i nterchangeabl y with horizontal wavelength. Frequency 
in this case relates to the time series of the in situ data. 
The true period or true frequency of a wave will be dis- 
tinguished from the period or frequency of the in situ data 
by calling the former the 'temporal" period or "temporal" 
frequency. 

The best method to obtain temporal charac ter i st i cs of 
the disturbances is to utilize data measured by a ground 
f ac i 1 i ty , such as an incoherent scatter radar or Doppler 
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ionospheric sounder, during a pass of the satellite. Then, 
from the satellite and ground measured data one can deter- 
mine both spatial and temporal characteristics, check the 
dispersion relationship and determine the amount of energy 
and momentum that they transport. 

Unfortunately, uery few satellite passes over ground 
facilities have been reported where data are available from 
both for the same time, suitable for wave analysis. Even 
when ground stations are operating, they may not be in the 
proper mode for determining the desire -» parameters. Prear- 
ranged measurements may not be successful too because of a 
lack of interesting events. 


In the absence of such corroborat i ve experiments, one 
may try to take advantage of the simultaneous measurements 
of many parameters by AE-C instrumentation to extract the 
temporal characteristics. This would be possible if the 
measured phase and amplitude relationships between the 
responses of individual parameters could be modeled. The 
simplest models are plane wave models. However, more com- 
plex modeling may be necessary; for example, a spherical 
model incorporating a dipole magnetic field. Here, we have 
attempted the simplest, Just plane wave modeling, and the 
results are reported. 


Satellite data for a full orbit were used. The data 
are typical of many orbits, and the results of the study are 
typical of results in analyzing other orbits of AE-C and 
also AE-E. The particular data set shown here was chosen to 
exemplify the results. The orbit is 5133 on January 20, 
1975. The magnetic field index k during the day w-s no 
greater than 3+ . Since the in ' i tu data were only for one 
altitude, vertical profiles can only be surmised by using 
modeling techniques. The data were available as 15 second 
samples. They come from the so-called AE-C 'Unified 
Abstract' <UA> files. 


The remainder of the body of this paper is as follows. 
Section 2. describes the data that were analyzed, the satel- 
lite orbit and geophysical conditions. Section 3. describes 
the techniques of data processing that were used to study 
the data. Section 4. contains results from studies of some 
statistical properties, such as auto- and cross correlation. 
Section 5. presents an attempt to model the responses of the 
geophysical parameters that were measured. Finally, the 
results of the study are discussed in Section 6. 


2. Pita frnt.lmd 
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Figure 1 is * graph of the densities o< atomic oxygen, 
molecular nitrogen, halium and tha electrons vs tima in 
seconds UT -for AE-C orbit 5133. Tha ordinate seal a is loga- 
rithmic. Density is in partielas par cm^ . Tha highast curve 
is for oxygen, followed in dacraasing dansi ty by nitrogen, 
halium and tha alactrons. Abscissa scalas also includa 
altitude, latitude, longitude, local tima, and tha magnatic 
field L shall paramatar <McI Iwai n , 1 961 ) . Tha altituda of 
tha orbit is in tha ranga 260 km ± 13 km. 15 sacond sam- 
plas ara plottad in tha figura. Thay ara so close that a 
continuous lina appears in almost all cases; parts of tha 
electron density curva show individual points where changes 
ara large. Tha nitrogan and halium curves ara coded with 
small crosses, so that thay ara heavy in appearance. Each 
curva contains 360 samples. One may note small intervals 
where there ara data drop-outs in tha individual curves. 
These occur simultaneously for tha neutral gas densities, 
since these parameters are measured by the same instrument. 
Fluctuations are apparent in all these densities, though 
soma of tha variations, in tha electron density are quite 
large and not necessarily connected with wave phenomena. 
For example, cn clrse examination one can see the so-called 
equatorial anomaly (Hanson and Moffett, 1966; Baxter and 
Kendal 1 , 1968) , as wall as rapid variations in tha auroral 
regions. 

Figura 2 is a graph of tha electron and ion tempera- 
tures for tha same period of time as in Figure 1. Altituda, 
longitude, L shall radius and geomagnetic inclination angle, 
in addition to time, ara shown as abscissa scales. Fluctua- 
tions ir. both temperatures ara quite evident, though those 
of the electron temperature are considerably greater. Tha 
electron temperature is more than twice tha ion temperature 
over most of the time span. There is a small interval of 
time at about 32298.8 seconds UT where the ion temperature 
appears double valued. A fluctuation between two values is 
actually occur ing hare, but the back and forth variation is 
obliterated by the ♦ symbol used to coda the trace. The ion 
temperature is higher in the southern hemisphere, which is 
for tha interval from 30952.5 seconds UT to 32298.8 seconds 
UT, than it is in the northern hemisphere for tha interval 
from 28260.0 seconds UT to about 30700.0 seconds UT. The 
large rise in the electron temperature within the first time 
interval is associated with tha southern auroral zona and 
appears broader and more intense than for the corraspondi ng 
parts of tha trajectory in tha northern hemisphere. 

A Mercator plot of tha orbit is shown in Figura 3. Tha 
times in seconds UT given on the abscissa scale of Figura 1 
ire shown in tha plot at about tha position of the satellite 
in its orbit *t these times. Tha direction of tha satellite 
in its orbit is shown by tha arrows. The data starts at 
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29260.0 seconds UT in the center with the satellite heaaed 
north over the eastern coast of the Un i ted State* on the 
right *ide of the map. The orbit crosses over the northern 
auroral region to Europe and A*ia reaching tfte right rage of 
the map at 40° ea*t longitude and about 40° north latitude. 
The trajectory transfer* over to the extreme left edge of 
the map where it i* shown heading south. In the southern 
hemisphere it is over the Indian and Pacific Oceans and 
passes close to Antartica in the southern auroral zone. The 
path turns north and the data stream and orbit both end in 
the center of the map just west of Mexico at 33645 seconds 
UT. 


A polar plot of the orbit in geomagnetic coordinates is 
shown in Figure 4. The path is shown as a thin line in the 
northern hemisphere and heavy in the southern hemisphere. 
The pole is the magnetic pole of the equivalent dipole 
field. The path shown is whet would be observed looking 
down over this pole if the earth were transparent. The 
position at selected times in seconds UT in the orbit are 
shown and may be related to various features evident in Fig- 
ures 1 and 2. Local times may be obtained from Figure 1. 

Magnetic and solar activity at the time were rela- 
tively quiet. The 3 hour planetary indices <Akasofu and 

Chapman, 1972) were *p*7,k«2, and the AE index (Davis and 
Suguira, 1966) was il4. p The solar 10.7 cm flux, F in W as 
77. 10,7 

3. Data Processi no 

Only the densities of oxygen, nitrogen, and the elec- 
trons and the ion temperature were processed to study their 
fluctuations, spectra and statistical correlations. This 
group of geophysical quantities will be referred to in 
places as the 'parameters' . Much of the data for helium and 
the electron temperature were considered too noisy for 
further study; some of these data are useful, but for much 
smaller time intervals. Though fluctuations of the parame- 
ters studied are evident in Figures 1 and 7, any relation- 
ship between these fluctuations is not clear. They be com-* 
clearer, however, on separating the fluctuations from the 
data. This separation was done by filtering the data, as 
will be explained below. 

Prior to filtering, it was necessary to close up any 
data g os; they are barely evident in Figure 1. This was 
done by fitting with a cubic spline interpolation scheme 
(Forsythe et a!., 1977). Data samples are then available 
for every interval throughout the entire time record. 

3 . 1_ Oa ta Normal i zat i on 


It was also necessary to normalize all data to be able 
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to compare different geophysical parameters. This had to be 
done prior to -filtering. Normalizing could be done in two 
way*. In the First, the data are normalized to the average 
value o-f the entire run o-f data. In the second, the data 
are normalized to the slowly varying component obtained by 
passing the data through a suitable lowpass Filter. This 
tatter process introduces some contamination at lower Fre- 
quencies and spreads any spectral Feature at higher Frequen- 
cies. Otherwise, there does not appear to be any diFFerence 
between the spectra oF data when processed in these two 
ways. The First methoo of normalizing was chosen to avoid 
any complications, and is presented here. 


Plots oF normalized densities oF oxygen, nitrogen and 
electrons and oF the normalized ion temperature are shown in 
Figures 5 through 8, respec t i ve 1 y . These are linear plots 
vs time in seconds UT . The entire group of data For each 
parameter are shown. Since each instrument was turned on at 
slightly difFerent times, there are small dIFFerences in the 
initial times in Figures 7 and 8 From those of Figures 5 and 
6 . All the records shown end at the same time. The fluc- 
tuations are far more evident in these figures than in Fig- 
ures 1 and 2. One can even see sane similar features in the 
curves of o«ygen and nitrogen, which become more apparent on 
superposing their graphs. Similar associations seem 
unlikely for the graphs oF the electron density and ion tem- 
perature . 

2.* 2 F i 1 ter i no 

The Fluctuations are extracted Fran the normalized data 
by passing the latter through a highpass Digital filter. 
The Filters employed are based on a design given by Hamming 
(1977). It is a nor.recurs i ve .'ilter which inc’udes a Lanc- 
zos window to reduce ripple* and minimize the Gibb's 
phenomenon. The Filter is described by its order and cutoff 
Frequency. Figure 9 is the transfer function of a typical 
Filter. Its order and cutoff frequency are 40 and 0.003 Hz, 
respec t i ve 1 y . The frequency scale in the F i gure is loga- 
rithmic and normalized to the sampling time, IS seconds, so 
that 0.003Hz corresponds to 0.045 on that s-.ale. The value 
of the Nyquist frequency on the same scale is 0.5, the last 
point on the abscissa axis. The value of the transfer func- 
tion at the cutoff Frequency is 0.5. Since the order is the 
number of weighted in-sequence samples that are summed to 
obtain the filter output, this many data are lost at the 
beginning and end of the record. Therefore, it is desirable 
to use the lowest order consistent with a reasonable steep- 
ness at the cutoff transition. The selection .of a filter 
order is a compromise between the number of data points 
available at the filter output and the suppression of fre- 
quency components i i the cutoff region. Thus, an order of 
40 was found satisfactory based cn trial and error and 


judgment of the quality of the power spectra of the output 
of the filter. All filters referred to, hereafter* are of 
this order. The original 360 data samples are reduced to 
282 samples at the filter output, which is adequate to pro- 
vide 256 data samples, a power of two, as needed by the Fast 
Fourier Transform. 

The choice of cutoff frequency is bas?d on the emphasis 
given to the various frequency components in the signal. 
For too low a cutoff frequency (0.001Hz and less), the tran- 
sition at cutoff is inadequate to keep out strong low fre- 
quency components without significant increases in the 
order. Furthermore, the filter cutoff must be selected so 
as to eliminate the fundamental frequency which has a period 
equal to the entire time interval of the data record, and, 
for safety, at least those frequencies up to the fifth har- 
monic of this fundamental. This part of the spectrum con- 
tains much power, but is of questionable validity because of 
the finite size of the data record. If not filtered, it 
would impair processing of the higher frequencies components 
of the spectrum. Because all data were found to decrease in 
almost a power law fashion with frequency, the lowest fre- 
quencies in the pass band are emphasized. Several different 
cutoff frequencies were used to take advantage of this 
feature, as will be explained. 

There was concern that spurious frequencies may be 
introduced by the Gibb's phenomenon as a result of the 
filter cutoff. Figure 10 demonstrates that such effects are 
insignificant. Three curves are shown in the figure plotted 
against time. The upper curve is the same as the nitrogen 
density curve in Figure 1. The middle curve is the output 
of a highpass filter with a cutoff frequency of 0.001Hz. 
The bottom curve is the same, but for a cutoff frequency of 
0.003 Hz. The curves have been lined up so that their 
abscissas in seconds UT match. Note the very excellent 
alignment of per turbat i ons in the original data with the 
fluctuations in the outputs of the two filters. No spurious 
frequency components are detectable in the outputs of the 
f i 1 ters . 

3«3 FI uc tuat i ons 

The fluctuations may be clearly observed at the output 
of a high pass filter. Figures 11 through 13 are plots of 
fluctuations vs time for a high pass filter of 0.003 Hz 
cutoff frequency. 0.003 Hz corresponds to a period of 333.3 
seconds, or a horizontal wavelength of 2600 km. The fluc- 
tuation data used to plot Figures 11 to 13 are normalized to 
their maximum values, so they that may be graphed between 
ordinates of -1.0 to 1.0. Each of the figures contains 
plots of the fluctuations of two parameters to study their 
correlation and relative phase. The fluctuation of oxygen, 
the major constituent, is used as the standard for 
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comparison in all three figures. Figure 11 compares oxygen 
and nitrogen density fluctuations. Figure 12 compares oxy- 
gen density and ion temperature fluctuations, and Figure 13 
those of oxygen and electron densities. In all three fig- 
ures the fluctuation of oxygen is shown as a solid line, and 
the other parameter as a dotted line. In Figure 13 the do*s 
are so c i ose in many places as to make the line solid. 
Nevertheless, the trace can be discerned. Oscillations are 
quite clear in all three figures. Dominant periods are evi- 
dent, adding construct i ve 1 y in some places and destruc t i ve 1 y 
elsewheres. It is also quite clear in all three figures 
that there are correlations between the fluctuations. That 
between nitrogen and oxygen densities is most outstanding. 
Oxygen density and the ion temperature in Figure 12 appear 
anti correl ated. The electron density and oxygen fluctua- 
tions in Figure 13 are certainly correlated with respect to 
the time durations of the individual oscillations. 

However, the electron density fluctuations are at first 
close to being in phase with those of oxygen, but from about 
30585.0 seconds UT to about 32235.0 seconds UT they are 
nearly out of phase. The electron density fluctuations 
appear relatively small in most places in this time inter- 
val, when normalized to the scale -1,1. The region of the 
equatorial anomaly is centered just about 30585.0, so that 
the satellite is passing from the northern to the southern 
hemisphere here. It seems that the electron density fluc- 
tuations from this filter, 0.003 Hz, do not relate to oxygen 
in the same way in both hemispheres. 

It is clear from the figures that the periods of the 
oscillations are about the same for all the parameters, even 
the electron density, and that more than one periodicity may 
be present. Measuring the spacing between several peaks for 
all parameters, one finds periods of the order of 140 
seconds and 250 seconds, correspond i ng to about 1090 km and 
1950 km, respectively, for a satellite velocity of 7.8 
km/sec. Ti.r shorter period wave is outstanding at time 
29760 seconds UT, whereas the longer period is readily 
observable in the fluctuations at about 30585 seconds UT and 
later. One may still observe the shorter period in features 
throughout . 

Maximum and root mean square values of the fluctua- 
tions are listed in Table I for oxygen, nitrogen, ion tem- 
perature and electron density for the 0.003 Hz filter. The 
rms fluctuations of nitrogen is about twice that of oxygen, 
whereas the rms fluctuations of the ion temperature is 
slightly greater than that of oxygen. The rms and peak fluc- 
tuations of the electron density are more than ten times 
that of oxygen, even though they appear small in Figure 13 
due to their being normalized to a scale of -1 to 1. 
Whereas the rms electron density fluctuations relative to 
the background density is about 20%, that of the neutral 


TABLE I 


Peak and Root Mean Square Fluctuations 
Filter Order = 40. Filter Cutoff Frequency = 0.003 Hz 


Ratio 


Parameter 

Peak Value 

RMS Value 

Peak /RMS 

Oxygen Density 

; 0525 

.0159 

3.3 

Nitrogen Density 

. 107 

.0338 

3.2 

Ion Temperature 

. .0728 

.0182 

4.0 

Electron Density 

1.255 

.203 

6.2 


TABLE II 



Peak and Root Mean Square Fluctuations 


Filter Order = 

40, Filter Cutoff 

Frequency = 0. 

006 Hz 

Parameter 

Peak Value 

RMS Value 

Ratio 

Peak /RMS 

Oxygen Density 

.0428 

.00926 

4.6 

Nitrogen Density 

.0756 

.0194 

3.9 

Ion Temperature 

.0809 

.0150 

5.4 

Electron Density 

.468 

.0867 

5.4 


11 


ORIGINAL PAGE IS 
OF POOR QUALITY 

densities and ion temperature are only a -few percent at 
most. Such large relative deviations of the electron der- 
si ty are well known from observations of TID's (Thome, 
1968). The ratio of peak to rms values range from 3 to 6 
for the four parameters. 

The location of the satellite in geomagnetic coordi- 
nates is also of interest relative to the graphs of Figures 
11 to 13. The key times in the plots are shown in the tra- 
jectory in Figure 4. It can be seen that the large oscilla- 
tions of nitrogen and oxygen in Figure 11 start when the 
satellite is close to the probable position of the auroral 
oval and continue to and past the geomagnetic equator at 
3058S seconds UT into the southern hemisphere. Oscillations 
are smaller in the southern mid-latitude zone, but pick up 
again as the satellite enters the region of the southern 
auroral zone between 31410 seconds UT and 32235 seconds UT . 
The ion temperature in Figure 12 behaves about the same way, 
but appears to have the strongest fluctuations in the 
auroral zone. The electron density fluctuations in Figure 
13 appear fairly uniform starting from 29760 seconds UT 
across the northern and- southern regions. There are some 
strong peaks on the northern side of the geomagnetic equa- 
tor, which is associated with the equatorial anomaly that 
may be clearly seen as the first of the two outstanding 
peaks near the middle of Figure 7. The southern part of the 
equatorial anomaly, the second and smaller peak in Figure 7, 
does not produce as large an excursion in Figure 13. It 
appears that the oscillations of the fluctuations obtained 
with the 0.003 Hz filter are associated more with northern 
hemispheric activity than is the case in the southern hemi- 
sphere, except for the hotter southern auroral zone (see 
Fi gure 2) . 

Figures 14 to 16 are the same as Figures 11 to 13, 
except that the filter cutoff frequency is increased to 
0.006 Hz (167 seconds period or 1300 km horizontal 
wavelength). Correlation in the fluctuations in all three 
figures is quite apparent. Even the electron density fluc- 
tuations in Figure 16 correlate better throughout, in con- 
trast with Figure 13; much more of the record is close to 
being in phase with the oxygen fluctuations, though phase 
shifts are a little larger in the equatorial region and the 
southern hemisphere. Nitrogen and oxygen appear to be in 
phase in Figure 14, and the i cn temperature in Figure 15 is 
shifted somewhat away from ant i -corre 1 at i on . In all figures 
there is a dominant fluctuation. In contrast with Figures 
11 to 13, the dominant oscillation in Figures 14 to 16 per- 
sists more clearly throughout much of the time span and 
appears to be on the order of the 140 seconds oscillation 
noted for Figures 11 to 13. This persistence is quite 
remarkable. Variations in amplitude also appear which imply 
the presence of waves of other periods, or horizontal 
wavelengths, that are beating against each other. 
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P»*k and rms valuta -for the fluctuations of Figures 14 
to 16 are given in Table II. Note that the values are 
smaller than those of Table I. The rms value's for nitrogen 
and oxygen are small er by about the same amount. The ion 
temperature fluctuations have not changed very much. The 
rms electron density fluctuations have decreased more than 
any of the other parameters. These changes result from the 
nature of the spectra which is found to decrease with 
wavelength; that is, there is more energy in the longer 
wavelengths, as will be shown in the next section. 

Figure 14 to 16 may also be analyzed with regard’to the 
position of the satellite. Much the same can be stated with 
respect to nitrogen and oxygen in Figure 14 as was stated 
for Figure 11. The comments relative to ion temperaature 
fluctuations in Figure 12 still pertain with regard to Fig- 
ure 15, and the same is true for the electron density fluc- 
tuations in Figures 13 and 16. Activity in the two hemi- 
spheres for scale sizes as indicated in Figures 14 to 16 
appear to be the same as for the larger scale sizes of Fig- 
ures 11 to 13. 


3.4 Spec tr al Anal ys i s 

Power spectra of the fluctuations are estimated using 
the Maximum Entropy Method <MEM) <Ulryck and Bishop, 1975) 
because of the finite size of the data sample and because of 
the lack of knowledge of data outside the time interval. 
The Fast Fourier Transform <fft) is also employed to help in 
selecting the order of the MEM and to judge its spectrum. 
Comparison with the periods of the fluctuations is neces- 
sary, as well, to insure that the MEM order and spectrum are 
consistent. 

Attempts were made to use the Akaike, or final predic- 
tion error <fpe), criterion to select the MEM order <Ulryck 
and Bishop, 1975). Figure 17 is a plot of the fpe vs. the 
MEM order for nitrogen fluctuations for a 0.0015 Hz filter. 
One chooses as the order the value that minimizes the fpe. 
Unf or tunatel y , the minimum is too broad and shallow to be 
useful for this purpose. The fpe is found to be similar for 
all other cases, so that the selection of MEM order is based 
on comparison with the periods of the fluctuations and with 
fft spectra. 

All spectra and fluctuations shown in this section are 
for nitrogen, and they are obtained with various filters. 
The spectra for oxygen density, ion temperature and electron 
density fluctuations are similar and need not be shown. In 
all cases a set of 256 fluctuation data are processed, all 
for the same time interval. 

An example of fluctuations is shown in Figure 18 which 
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is for a high pass filter with a cutoff of 0.0015 Hz. The 
ordinate is scaled for the actual magnitude of the fluctua- 
tions in contrast with the plots of Figure 11 through 16. 
The periods, though not exactly constant, are quite evident 
and are found to be about 500 seconds for one set of fluc- 
tuations and 250 seconds for another. These times 
correspond to horizontal wavelengths of about 3900 km and 
1950 km, respectively. Figure 19 is a plot of the fft spec- 
trum of the fluctuations in Figure 18. The abscissa scales 
are the frequency in millihertz (mHz) and the wavelength in 
km. The frequency scale extends from 0 to the Nyquist fre- 
quency for 15 seconds data samples, 0.0333 Hz. Spectral 
peaks are labeled by numbers that represent harmonics of the 
fundamental frequency which has a period equal to the total 
time span of the 256 samples in Figure 18, namely, 3825 
seconds. These peaks are at harmonic numbers 4, 8, 10 and 
16. The peak at harmonic number 4 is ignored as too low, 
since it may be affected by the finite data sample. The 
peaks at harmonics 8 and 16 correspond to periods of 478 and 
239 seconds, or wavelengths of 3729 and 1865 km, respec- 
tively. These waves confirm the periods evident in Figure 
18. The peak at harmonic 10 may be real , or i t may be the 
result of the variations in the 500 second period in Figure 
18. Since the computed spectrum is at discrete frequencies, 
the broadening associated with this variation may appear as 
a nearby additional, but false peak. 

One must select an order for the MEM. If too low, sig- 
(ii 4 leant spectral features may be submerged. If too great, 
spurious features may arise. It is difficult to judge what 
order to use. The fpe in Figure 17 is not very helpful. A 
value anywhere from 4 to 16 appears to produce a shallow fpe 
minimum. For an order of 15 (not shown here) the MEM spec- 
trum misses the waves of harmonic number 8 and 16 entirely. 
Only one peak arises at harmonic 11. Increasing the order to 
20 results in a spectrum with the desired peaks. This spec- 
trum is shown in Figure 20. The two largest are for harmon- 
ics 8 and 16. Another much smaller peak appears at harmonic 
numbers 30, and still smaller peaks are suggested at har- 
monic numbers 42 and 60. A logarithmic plot of the spectrum 
is shown in Figure 21 which is for frequencies from the 3th 
harmonic to the Nyquist frequency. The third abscissa scale 
in the figure gives the harmonic numbers. The peaks at the 
higher frequencies, or shorter wavelengths, can be seen in 
this plot. Note the power law fall off with frequency. 
The straight line is a least square fit to the spectrum. 
The slope of the line on a log-log plot is -3.1, which means 
a power law with frequency varying as f~3«l over at least a 
16,5 :1 frequency range. 

By passing the normalized data through high pass 
filters of higher cutoff frequencies which filter out the 
lower, more powerful parts of the spectrum, one can get a 
better sense of the higher frequencies or shorter wavelength 


components in the data. Figure 22 is a plot of the MEM spec- 
trum of order equal to 20 -for nitrogen using a -filter with a 
cuto-f-f of 0.003 Hz. There is one dominant peak at about har- 
monics 15 to 16. The cutoff is sufficient to eliminate the 
peak at harmonic 8 observed in Figure 20, since 0.003 Hz 
corresponds to a non-integral harmonic of 11.5. This cutoff 
is too low to highlight waves of wavelength shorter than 
that of harmonic 16. A filter with a cutoff of 0.006 Hz will 
do this by eliminating the peak at harmonic 16. Figure 23 
is a plot of the fluctuations at the output of this filter. 
Careful measurements of the periods between peaks there give 
periods of about 170 seconds, 127 seconds and 106 seconds 
with about a ± 5% error. These correspond to harmonic 
numbers 23, 30, and 36 with harmonic number errors of ±1 at 
23 to ±2 at 36. Only one of these components are evident in 
Figure 21, that for harmonic 30. This component is the more 
persistent periodicity in Figure 23. The correspond i ng MEM 
spectrum of order 20 of the fluctuations in Figure 23 is 
shown in Figure 24, where we observe a strong peak at har- 
monic 31 and a weaker peak at harmonic 25. The latter con- 
firms the longer period in Figure 23. The period for har- 
monic number 36 is not seen in this spectrum. On increasing 
the order of the MEM to 35 with a filter cutoff of 0.0015, 
one obtains the logarithmic spectrum shown in Figure 25. 
Here we find peaks close to all expected periods based on 
analysis of Figures 18 and 23. The corresponding fft spec- 
trum (not shown here) has peaks at the desired values as 
well. 


On increasing the cutoff frequency further, we obtain 
fluctuations for nitrogen shown in Figure 26 which is for a 
filter cutoff of 0.0095 Hz. Careful scaling, with account 
taken of measurement accuracy, suggests periods for harmonic 
numbers 34 to 38, 39 to 43, and 57 to 63. Peaks at, within, 
or close to these values are exhibited in the spectrum of 
Figure 25. The fft spectrum (not shown here) has many more 
peaks for this case, some of which may be extraneous. Fluc- 
tuations for a filter with a cutoff of 0.013 is shown in 
Figure 27. The periodic wave there, with error accounted 
for, is in the harmonic range 60 to 66. This wave may 
account for the peak at harmonic number 63 in Figure 25. The 
fft spectrum for this case, not shown here, exhibits many 
peaks that must be extraneous. 

Since most of the waves suggested by the nitrogen den- 
sity fluctuations appear to be exhibited in the MEM spectrum 
for the order 35, this order and spectrum are taken as ade- 
quate for the purpose here. Similar tests with different 
filters, as used with nitrogen density data, were employed 
with the other parameters. These tests confirmed the use of 
an MEM order of 35 for these parameters. Thus, this order 
is used to obtain the MEM power spectra for the fluctuations 
of all of the parameters after passing their data through a 
high pass filter with a cutoff frequency of 0.0015 Hz. A 
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comparison against harmonic numbers of peak frequencies 
these parameters is given in Table Ilia. The 
wavelengths correspondi ng to the oxygen spectral 
also given in the table. The power 
least square straight line -fits to the 
are given in Table 1 1 1 b . 


■for 

hor i zontal 
peaks are 
law exponents -for the 
logarithmic spectra 


The MEM tends to broaden spectral peaks. Nevertheless, 
Table Ilia indicates spectral agreement -for as many as three 
or -four parameters -for a number o-f wavelengths, 1845 km, 
1193 km, 942 km, 523-533 km and 444 km. In addition, oxygen 
and nitrogen spectra correspond -for 3729 km. The i ora tem- 
perature may also have the. same -feature at this wavelength, 
though its value is somewhat o-f-f. The ion temperature and 
oxygen have -features at 728 km, but nitrogen and the elec- 
tron density spectra both peak at a smaller wavelength in 
the range 478-493 km. Similar -features also appear possible 
■for shorter wavelengths at 409 km and less. Still shorter 
wavelength spectra may be limited by noise. 


From the numbers in Table Ilia, it appears that the 
spectra of neutral densiti.es correspond for many wavelengths 
from the order of 4000 km to hundreds of kilometers. 
Correspondence of ionization with neutral spectral features 
over many thousands of kilometers may not really occur, 
since plasma is subject to many other physical processes, 
particular those involving electric fields that do not 
affect neutrals. Nevertheless, the spectral feature at 1845 
km is present in the electron density spectrum as well as 
those of the neutral densities and the ion temperature. 
Since the ion temperature is coupled to the neutral tempera- 
ture at an altitude of 240 km, it would not be surprising to 
find that the ion temperature spectrum has features at very 
long wavelengths like those of the neutral gases. 

The power law type of decrease with spatial frequency, 
or with decreasing wavelength, is present in all parameters. 
The exponents are different, as shown in Table Illb, though 
those of oxygen and nitrogen are close in value. The ion 
temperature spectral decay rate is smallest, whereas that 
for the electron density is faster than oxygen and nitrogen. 
Such power law decreases may arise from the sources as well 
as the medium, but may also be related to turbulence gen- 
erating processes. The reason for the different exponents 
is far from clear and requires further study. These charac- 
teristics are in agreement with ideas of Yan 2andt <1982). 


£• Statistical Aspects 

Statistical investigations here include only computa- 
tion of au tocorre 1 at i on and cross correlation functions and 
an au toregress i ve <AR) analysis (Jenkins and Watts, 1948). 
These quantities and functions require stationary statis- 
tics. Furthermore, the MEM technique requires at least 


TABLE III 


Possible Spectral Harmonics for 256 Data Samples 


Oxygen 

Density 

Nitrogen 

Density 

Ion 

Temperature 

Electron 
- Density 

Horizontal 
Wavelength (km) 




6 


8 

8 

7 


3729 

16 

16 

16 

16 

1865 

25 

24-25 

26-27 

24 

1193 

31 

31 

31 

32-33 

962 

m 

37 

- 


- 

41 

43-44 

41 

43-44 

728 

56-57 

56 

56 

55 

523-533 

64 

63 

64 

63-64 

466 

73 

73 

71 

71 

409 

82 

83 

81 


364 

90 

90 

92 

88 

332 


Power Law Exponents 
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quasi -stat i onar i ty to be meaningful. 

Sine* the number o-f data samples are limited, tests for 
stationarity were not attempted. A qualitative evaluation 
was used instead. For one, fluctuation signals, as in Figure 
11, were visually examined to determine whether the nature 
of the signal was maintained throughout the record. The per- 
sistent character of the records in Figures 11 to 16 suggest 
some degree of stat i onar i ty . 

Typical ensemble samples were computed using an autore- 
gressive analysis to check further for stat i onar i ty . If a 
sample has character i st i cs similar to those of the actual 
fluctuation signal, the autoregress i ve process was judged to 
be satisfactory and an indication of stat i onar i ty . This 
indeed was found to be the case for the four parameters that 
were studied. Other data not shown here for helium, a 
parameter that was not analyzed, were deemed quite noisy and 
poor for analysis. AR models for these data were vastly dif- 
ferent from the actual fluctuations, so that the helium data 
were judged unsuitable for a stationary process. 

Figure 11 exhibits both the nitrogen and oxygen flucta- 
tion signals. Figures 28 and 2 9 are AR models of order 20 
for oxygen and nitrogen, respect i vel y . Autocorrelation coef- 
ficients needed to compute the model are based on data 
passed through a 0.003 Hz filter. Note the similarity in 
the oscillatory character of the models to the actual fluc- 
tuations in Figure 11. The spacing between peaks in tnat 
figure appear narrower than those of Figures 28 and 29, but 
this is due to the smaller time interval in the abscissa 
scale in those figures. Figure 11 is 300 seconds longer in 
time than Figures 28 and 29. The AR model for each neutral 
gas is but one ensemble samp '• e . Since the process is con- 
trolled by a noise source, a whole ensemble class may be 
generated. It is found that the characteristics persist, 
from which we conclude that the data are at least quasi- 
stationary. The same is found for the ion temperature and 
the electron density. 

Figures 30 and 31 are the autocorrel at i on functions of 
the oxygen and nitrogen fluctuations shown in Figure 11, 
respect i vel y . The oscillatory nature of these graphs are 
quite evident. The plots are against time lag intervals, 
where one lag represents 15 seconds. The oscillations per- 
sist out to at least 80 lags in both cases. Further, the 
oscillations in these two graphs are quite similar. The 
periodicity is almost the same. Also, other frequencies 
must be present to account for the nature of the variation 
with the time lag. The separation between the peak at zero 
lag and the next peak is scaled to be 267 seconds and 247 
seconds for oxygen and nitrogen, respec t i ve 1 y . These times 
agree with the periods in Figure 11 and the spectral peak in 
Figure 22. Small differences may arise from beats due to 
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The cross correlation function (CCF) between the fluc- 
tuations of two parameters is useful to obtain a sense of 
the degree of correlation and the phase shift between the 
parameters; that is, the phase angle of the transfer func- 
tion between the responses for the latter. Oxygen density 
fluctuations are used for this purpose to compare with the 
fluctuations of the other three parameters obtained using 
the same filter. Oxygen is chosen for comparison, since it 
is the main neutral constituent. Figures 32 and 33 are cross 
correlation functions plotted us. 15 second lags for*oxygen 
and nitrogen density fluctuations obtained with filters of 
cutoff frequencies of 0.003 Hz and 0.006 Hz, respect i uel y. 
Note the oscillatory behavior in both figures. Note too that 
the peak values of the CCF is quite adequate for correlation 
in both cases, being just more than 0.7 for the 0.003Hz 
filter and about 0.69 for the 0.006 Hz filter. The phase 
shift is nearly zero; there appears to be a slight phase 
angle with nitrogen leading (to the left of zero lag) for 
both filters. The best that can be determined is that the 
lead is about 1/2 a lag interval, about 7.5 seconds. The 
oscillatory behavior persists on both sides of zero to 
longer lags than shown. The spacing between peaks in Figure 
32 on either side of the large peak near zero lag is 16 1/2 
lags. Thus, the dominant period is 247.5 seconds which 
corresponds to an on-track horizontal wavelength of 1931 km 
for the fluctuations from the 0.003 Hz filter. This period 
agrees with the correspondi ng values from the autocorrela- 
tion functions in Figures 30 and 31, the spectrum in Figure 
22 and the estimated period in Figure 11. The dominant 
period for the 0.006 Hz filter in Figure 33 is 8 1/2 lags, 
or 127.5 seconds, which corresponds to a horizontal 
wavelength of 995 km. This period and wavelength are 
roughly in accord with the spectrum in Figure 24 for nitro- 
gen, the harmonic 31 listed in Table Ilia and the periods of 
the fluctuations in Figure 14. The 1/2 lag interval phase 
lead corresponds to a leading phase angle of 11° for the 
0.003 Hz filter and 21° for the 0.006 HZ filter. The peak 
CCF becomes even larger for fluctuations of nitrogen and 
oxygen from filters of cutoff frequencies as high as 0.013 
Hz. Oscillatory behavior of the CCF persists with periods 
that are shorter than for filters of lower cutoff frequen- 
cies. 


The CCF for oxygen density and ion temperature fluctua- 
tions show some correlation but not as strong as that 
between nitrogen and oxygen. That for a 0.003 Hz filter is 
shown in Figure 34. This CCF is quite irregular in contrast 
with that for Figures 32 and 33. Note that the peak is nega- 
tive with a value of -0.38, implying an t i -corre 1 at i on . The 
period cannot be determined from Figure 34. Figure 35 is a 
plot of the CCF for oxygen and ion temperature fluctuations 
from a 0.006 Hz filter. There are more oscillations than in 
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Figure 34 ; part i cul arl y on the leading aide, and the peak is 
nearly ant i -corre l ated w « t to a value of -0.41. The negative 
peak leads by somewhat less than 1 lag interval. The plot 
suggests the presence of more than one wave with resultant 
interference. Such interference would produce an error in 
scaling the dominant period, as well as the phase shift. 
The spacing of peaks on the lead side is 10 lag intervals, 
whereas that on the other side is 7 lag intervals. The mear 
of 9 1/2, if meaningful, agrees with that for oxygen and 
nitrogen in Figure 33. The maximum negative peak at -1 lag 
on the leading side in Figure 35 corresponds to a lagging 
phase angle of -138° for a period of 8 1/2 lags, whersas the 
phase angle would be -159° for a maximum negative peak at a 
-1/2 lag interval. Thus, the phase shift appears to lie 
between the two values, -138° and -159°. The CCF peak 
reaches -0.49 for a filter cutoff of 0.008 Hz (not shown 
here), where an t i -corre 1 at i on persists. The cutoff for this 
last filter corresponds to what would be a non-integral har- 
monic value of 30.6, or a horizontal wavelength of 975 km. 
Thereafter, at still higher filter cutoffs the CCF between 
oxygen and ion tempe ature becomes smaller, but the correla- 
tion to nitrogen remains stronger, reaching -0.46 at a 
filter cutoff of 0.010 Hz. This suggests the possibioity 
that, for the shorter wavelengths, the ion temperature is 
coupled to the nitrogen density more than the oxygen den- 
sity. 


Though the degree of cross correlation of ion tempera- 
ture and oxygen is not large, ant i -correl at i on stands out, 
and for wavelengths of 1000 km and less a periodicity is 
evident in the cross correlation function. This periodicity 
relates to that between nitrogen and oxygen and suggests 
that the cross correlation of the ion temperature is signi- 
ficant though the peak values are less than 0.5. 

Figure 36 is a plot of the CCF of oxygen and electron 
density fluctuations for a 0.003 Hz filter. Though the peaks 
are not large, an oscillatory variation is quite evident. In 
addition, there appears to be interference which may be from 
another wave that makes it difficult to dete r mine the true 
peak near zero lag. The lag separation between peaks for 
negative lag intervals is about 17 whlcn is close to the 16 
1/2 found for nitrogen and oxygen for the sane filter. The 
separation between the other peaks is much smaller. The 
changes in the phase shift of the fluctuations, as described 
in Section 3.3 relative to Figure 13, may be the reason for 
these difficulties and the poor magnitude of the *<>ak CCF. 

Figure 37 is the CCF of the electron density and oxygen 
for a filter with a cutoff of 0.006 Hz.- Much more osc i 1 I a- 
tory behavior is evident. The largest peak is well to the 
right of the peak near zero lag. This suggests interference 
between wave components. The peak at zero lag is 0.39, 
whereas the peak to its right is 0.45. The interval between 
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the peaks about zero lag is 9 1/2 
larger than that found in Figure 33 for nitrogen and oxygen. 
The discrepancy may arise from i n ter f erence Nevertheless, 
the cyclic variation suggests better cross correlation than 
indicated by the magnitudes of the peaks. 

Since neither the ion temperature nor the electron den- 
sity CCF has peak magnitudes of 0.5 or more, they would not 
be expected to correlate well with *«cn other, as their own 
CCF does indicate, though not shown here. 

There is a small phase lead for the CCF of Figure 37 of 
about 1/2 lag. This would correspond to a phase lead of the 
electron density fluctuations relative to the oxygen fluc- 
tuations of 21°, No reliable phase shift may be obtainec' 
from Figure 36 for the 0.003 Hz filter. As the filter fre- 
quency becomes larger, the cross correlation of electron 
density and oxygen density is found to become poorer, 
undoubtedly because of greater noise at the higher frequen- 
cies. Figure 37 suggests correlation between electron den- 
sity and oxygen for waves of the order of 1000 km. Correla- 
tion at still smaller wavelengths may be obscured by noise, 
though indicated in Table Ilia. 

One may use the crude phase shifts from the CCF calcu- 
lations and the mean square fluctuations of the parameters 
to describe transfer functions between the responses of the 
parameters for further use in modeling. This is described 
next . 
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3. Thtore t i cal Mode 1 i n q 

Tht presence of waves is clearly shown in Figures 11 to 
16. Mort than on# wave is prtstnt as shown by tht spectra 
of Figurts 20,21,22,24 and 26. Tht waves that dominate 
dtptnd on tht filttr, sinct tht sptctrum dtcrtasts with 
dtcrtasing wavelength as shown in Figurts 21 and 25 and as 
required by tht various sptctral txpontnts listtd in Tablt 
1 1 1 b . Tht iowtst frequency compontnts passtd by tht filttr 
dominatt. Two dominant ptaks may bt obstrvtd in Figurt 24 
for tht 0.006 Hz cutoff filttr. Tht first ptak is at a 
wavtlangth of ’193 km, whtrtas tht stcond is at .962 km. 
Thtst two ptaks would btat against tach othtr to product a 
btat patttrn. Ont may inttrprtt tht tnvtlopt variation in 
Figurt 23 as arising from such a btat patttrn. Of courst, 
wt must bt cartful in inttrprtting thtst fluctuations and 
sptctra, sinct wt art dtaling with a finite data samplt that 
is just a snapshot of tht nitrogtn distribution at ont alti- 
tudt in this cast. This distribution is tht rtsult of sto- 
chastic sourcts and a numbtr of geophysical processes that 
art at work and may not actually rtsult from the beating of 
two wavts. 

It is important to know whether the wavts art traveling 
or stationary, or both. If traveling, they art transporting 
energy and momentum from their sourcts. If standing wavts, 
they may bt stationary spherical modes due to ont sour e 
<Mayr tt al , 1982), or they may originate from interfere,' .# 
of waves from near duplicate stochastic sources in .j.'th 
auroral zones, such as could arise from diffuse aurora. 
This type of auroral activity may produce the same effects 
in both hemispheres (Lui et al , 1982; Akasofu, 1974). 
Standing wavts would die out by dissipation simultaneously 
over a large region, much like plucking a violin string. It 
has been noted in Section 3.3, for example, that there is 
more activity in the level of the fluctuations in the entire 
northern hemisphere portion of the trajectory than the 
southern hemisphere, with the exception of the southern 
auroral zone. It may be that these waves are primarily due 
to one source and both stationary and traveling waves are 
prtstnt. Analysis of simultaneous ground measurements of 
ionization and satellite measurements of in situ neutral and 
ionization data may aid in understanding better these wave 
systems. 

In the absence of such ground measurements, and with 
the limitations of in situ measurements, we can only test 
for traveling waves. To determine whether the waves are 
pure traveling wavts we assume a plane parallel model and a 
plane wave of some horizontal wavelength and temporal fre- 
quency. Tht assumption of a plane parallel model is a crude 
approximation for propagation over a sphere. It still may 
bt useful, sinct, as shown by Francis (1972), the waves tend 
to follow the curvature due to rotation of the direction of 
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gravitational acceleration as the wave propagates. One 
seeks the values oF wave parameters that provide the closest 
duplication of the ratios oF the responses «oF nitrogen den- 
sity, electron density and ion temperature Fluctuations to 
that oF oxygen; that is, the transFer Functions between 
these various geophysical parameters. In this way we can 
obtain the wave properties not directly observable From the 
spatial distribution, such as the temporal Frequency and 
direction ' propagation. The experimental transFer Func- 
tions are given by a magnitude and phase angle. The magni- 
tude is obtained From the ratio oF the root mean square 
Fluctuations oF the geophysical parameters (Figures 11 
through 16), and the phase angle is estimated From the cross 
correlation Function (Figures 32 through 37). Such transFer 
Functions ari determined For each Filter. We give here two 
examples, that For 0.003 Hz Filter and For the 0.006 Hz 
Filter, to be consistent with the previous material. Table 
IV lists these transFer Functions. 



* positive angles represent leading phase angles, whereas 
negative angles are lagging phase angles. 

♦ based on periodicity in cross correlation Function in 
Figures 32 and 33. 


There is, oF course, a question oF the meaning oF these 
transFer Functions, since the amplitudes oF the Fluctuations 
vary over the time span oF the data. One can see in Figures 
11 through 16 that, although the amplitude oF the Fluctua- 
tions vary and there are diFFerences in their variation, the 
Fluctuations oF all the parameters behave somewhat in the 
same way. ThereFore, the ratio oF the rms Fluctuations may 
be useFul within some limits For the purpose here. Simi- 
larly, the cross correlation Function may be an acceptable 
indicator oF the phase shiFt, since ther? is some conFirma- 
t i on with phase relationships observable in Fiyures 11 to 
16, though perhaps only partly so For the electron density. 
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In any event, adoption of these magnitudes and phase shifts 
may be the best one can do with the data. 

* 

The mode) employed assumes two species, one major in 
density, the other minor. Then all minor species do not 
seriously affect the major species, and one may have any 
number of minor species. Each is treated in turn as the 
second species together with the major constituent. It is 
assumed that a gravity wave is moving through the main atmo- 
sphere and drags each of the minor species by collisions. 
This is a reasonable model -for ionization and oxygen. The 
nitrogen density is smaller than the oxygen density (Figure 
1) throughout the entire data span. Though nitrogen is not 
necessarily truly minor throughout, the approximation is 
still useful considering the coarseness of the transfer 
functions in Table IV. 

The model must also incorporate the earth's magnetic 
field to determine the response of the ionization. It is 
assumed that charged particles are constrained to move only 
along the magnetic field. Figure 4 shows the plot of the 
orbit in geomagnetic coordinates. The times along the 
abscissas in Figures 11 to 16, 18,23,26 and 27 are shown 
there. A major portion of the orbit for the times in the 
above figures (not the full time span shown in Figure 4) is 
at geomagnetic latitudes of magnitude greater than 30°, so 
that the magnitude of the dip angle exceeds 49° .for most of 
the orbit. The average of the magnitude of the dip angle is 
60° over the entire time span. Since the Tiode 1 employs a 
plane stratified geometry, it is not possib) to vary the 
dip angle very readily as could be done i t h a spherical 
model. Therefore, constant dip angles are used. The 
transfer functions for a value of 60° were compared with 
those for smaller dip angles. Since it was found that vari- 
ations due to the dip angle over a range of values were 
mostly secondary, the choice of a single angle for a given 
set of calculations is useful. However, it makes a differ- 
ence whether propagation is towards the equator or away from 
the equator. The calculations for zero dip angle yields 
results that are like propagation away from the equator with 
non-zero dip angles. The closest approximation to the 
values in Table IV were found for propagation towards the 
equator. Results for electron density are presented for dip 
angles of 20° and 60°. 

Since collisions between species is of importance, one 
must have adequate esti ates *)f the collision frequencies. 
The value of the collision frequency between nitrogen and 
oxygen, 0.372 sec , had been obtained from the hard sphere 
model (Chapman and Cowling, 19 ) for a mean temperature of 
830°K (as determined from the ion temper ature ) , a total den- 
sity of 10 y particles/cm 3 , and collision diameters of 0.18 
nm and 0.19 nm for oxygen and nitrogen, respec t i ve 1 y . Then, 
computations were made with this collision frequency and 
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values differing by -factors o-f two above and below to com- 
pensate -for differences in temperature and inadequacies of 
the hard sphere model. Ion-neutral collision frequencies 
are based on an expression given by Stubbe <1968) for 0 + , 
the principal ion, and atomic oxygen. The average ion tem- 
perature, 830°K, and a neutral density of 10^ particles per 
cm 3 were employed and gave a nominal value of 1 .75sec -1 . 
Acain values differing by factors of two were employed to 
compensate for inexact temperatures and inadequacies in the 
collision frequency estimates. 


5.1 Formul at i on 


The theory of coupling of species in the atmosphere to 
gravity waves moving through the entire atmosphere has been 
developed for two species <Gross and Eun, 1976 and 1978; Del 
Genio et al,1978 and 1979), though Mayr et al <1978) has 
modeled the problem for a number of constituents on the com- 
puter. In Gross and Eun , the second species is taken as 
minor relative to the first, with the advantages as 
explained above for multiple minor species. The theory of 
Del Genio et al is more general; it is also for two consti- 
tuents, but they may be of equal importance.. The formula- 
tion for nitrogen and oxygen densities used here are based 
on the theory of Gross and Eun. Neither of the above refer- 
ences treat ionization in an oxygen dominated atmosphere. 
There have been treatments of ionization density based on 
the continuity and momentum equations <Hooke, 1968), but 
none incorporate the ion and electron temperatures as fluc- 
tuating parameters. A formulation for ionization driven by 
a gravity wave in the neutral medium, with a superposed mag- 
netic field and i ncorporat i ng the ion and electron tempera- 
ture fluctuations, has been developed and used here. The 
details will not be given here for brevity. Only the 
expression for the transfer function will be given. Though 
the theory is not restricted to isothermal conditions, this 
assumption was made for simplicity. The vertical tempera- 
ture gradient is not known, and probably not important here. 
It would represent a refinement in the calculations. The 
formulation also neglects viscosity, but this may be impor- 
tant. It also ignores possible reflections from below 
<Franc i s, 1974) . 

The ratio, T^, of the nitrogen density r' ponse to that 
of the oxygen density, for a plane wave exp<-.a)t+ik x+ik z), 
may be shown to be given by equation <l)s x 2 
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M=ratio of the particle masses of n i trogeTf" '^b 1 
«28/16>j w is the temporal frequency, W, relative 
buoyancy frequency N, that is, w=u/N; K x 4ri( j k are 

numbers k^ and k_ in the horizontal and vertical 
respectively, normalized to the scale height of the 


Here , 

O 


oxygen 
to the 

the wave 

di rect i ons, 
main 


neutral constituent, H , so that K x =k x H 0 and K,=k z H iY is 
the ratio of specific heats of atomic oxygen <=5/3> ; w” n j s 
the col i si on frequency between oxygen and nitrogen, 
relative to the buoyancy frequency. The model assumes an 
isothermal atmosphere, so that the buoyancy frequency is 
given by 


N=''4F-7a/a 
~ n 


< 2 ) 


where g is the acceleration of gravity <?.05m/sec at 260 km 
altitude) and a n j s the velocity of sound. If w and k x , or 
w an< * **x * are rea1 » k 2 , or K z , will be complex in general. 
K 2 is related to w and K x by the dispersion relationship 
for an isothermal atmosphere: 


* 

K 2 s*±1/w £i lw *-<K 2 *l/4)w 2 + K 2 J - i /Z (3) 


The plus sign is for upward phase propagation and the minus 
sign for downward propagation for time harmonic variation 
exp<-iwt>. It is assumed in these calculations that the 
source is below the observation point, so that phase pro- 
pagates downwards, or energy propagates upwards. 

The ratio of the electron density fluctuations, R. , and 
the ion temperature fluctuations, T. , to that of oxygen is 
given by solving the matrix equation ^4) : 
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In <4> is the ratio of ion velocity along the field line 
to the acoustic velocity of the principal neutral consti- 
tuent; w, , K 2 and Y have been defined; C b and S b are 
short hand for cos 6 . and sin$ b , where 9 b is the dip angle, 
positive in the northern hemisphere and negative in the 
southern hemisphere; w. is the ion-neutral collision fre- 
quency, V. , relative to the buoyancy frequency, N. The 
ratio of electron to ion temperatures is a . It is assumed 
that the relative electron and ion temperature fluctuations 
are also in the same ratio. Electron temperature enters 
effectively only through the electron pressure, since colli- 
sions between electrons and ions are neglible at 2<£0 km 
altitude. Since coupling of neutrals to the electrons is 
primarily by means of the ions at this altitude, the assump- 
tion appears reasonable. The value of a tnat was chosen is 
based on the mean electron and ion temperatures (Figures 2). 
A value of 2 was adopted based on their averages. The quan- 
tity H. in equation (4) is the ratio of the ion density 
scale height to the scale height of oxygen. Assuming a beta 
Chapman layer at 260 km altitude, one can show that H.- 4 / 3 . 
(In this regard it should be noted that the satellite was in 
sunshine for over 707. of the data.) Other values of H. &r -e 
also possible, though it not expected that significantly 
different results would be obtained. 
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The coordinate system assumes 
northward, so that K j s positive -for northward propagation 
and negative -for southward propagation. 

5.2 Resu 1 ts of Mode 1 Cal cu 1 at i ons 

The ratio o-f the response o-f nitrogen to that -for oxy- 
gen using equation <1), and the ratios o-f the response of 
the electron density R. and ion temperature T- to that of 
oxygen, on inverting the matrix equation <4), were calcu- 
lated as a function of frequency for a given temperature T , 
n i trogen-oxygen collision frequency v t ion-neutral coll?- 
sion frequency Vj n , dip angle and horizontal wavelength 
X . These ratios will be referred to as transfer functions. 
Tne purpose of the calculations was to find those parameters 
which come closest to matching the values in Table IV. The 
normalized vertical wavenumber for each value of frequency 
is obtained from <3). Several values of both nitrogen- 
oxygen and ion-neutral collision frequencies were used. 
Values for the n i trogen-oxygen collision frequency were 
0.0?3sec , 0.186sec~ , 0.372sec~ 1 and 0.75sec“ . Values 

for the ion-neutral collision frequency were 0.44sec“*, 
0 ,87sec” , 1 ,75sec“* and S.Ssec"* 1 . Dip angles from 0° to 
±80° were employed. 

It was quickly found that waves traveling towards the 
poles did not provide suitable phase shifts or amplitudes as 
given in Table IV. Waves traveling towards the equator, 
however, did provide the phase shifts over certain ranges 
of the parameters. Emphasis in the calculations was there- 
fore placed on waves of this type, particularly for a dip 
angle of 60°, the mean magnitude. Temperatures other than 
830°K were also used; namely, 6Q0°K and 1200 C K. 

An example of the calculated nitrogen to oxygen 
transfer function amplitude and phase shift, as given by 
equation <1>, is shown in Figure 38 for a horizontal 
wavelength of 2000 km, collision frequency V =.75sec -1 and 
T q = 830°K. This plot is against the actual temporal fre- 
quency in mHz. The left scale is the phase shift as given 
by the dotted curve. The right scale is for the amplitude 
as given by the solid line. The points of the phase curve 
matching the phase shifts of the fourth column of Table IV 
are indicated by the 3's for the 0.003 Hz filter <ll°±i°) 
and by the 6 ' s for the 0.006 Hz filter <21°±1°). The ampli- 
tudes for this case exceed somewhat the values given in the 
third column of Table IV. Similar examples for the electron 
density and ion temperature transfer functions are shown in 
Figures 3? and 40 for a horizontal wavelength of 1000 km, 
ion-neutral collision frequency v . t g 7 sec - * and dip angle 
$l=60 o . The 3"s and 6's shown there are at the points of 
the spectra where the phase angles in columns 6 and 8 of 
Table IV are obtained. Since the phase shift for the elec- 
tron density fluctuations after passing through the 0.003 Hz 
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■filter is uncertain, no 3's are shown in Figure 39. A small 
phase shift is expected which would be smaller than the 
ualue indicated by the 6 ' s in the figure, so 4 that the fre- 
quency would be scrie ualue less than the frequency at the 
6"s. The amplitudes of the ion temperature transfer func- 
tion in Figure 40 are nearly unity in all cases, different 
from the values in column 5 of Table IV. The magnitudes of 
the electron density transfer function significantly exceeds 
unity. That shown n Figure 39 near the 6 ' s is about 17, a 
ualue greater than the values in column 7 of table IV. 

The resonance like peak in the electron density 
transfer function is due to a complex pole in the transfer 
function arising from coupling between gravity waves and a 
sort of ion diffusion wave. This type of resonance was 
described by Gross and Eun <1978) for helium. 

Figures 41 and 42 contain summaries of the calculations 
for a temperature of 830°K. These are typical of results 
for the other temperatures used in the calculations. Fre- 
quency in millihertz is plotted vs horizontal wavelength in 
each figure. The horizontal wavelength here is assumed to 
be the true wavelength, so that it takes values less than or 
equal to the maximum on-track wavelength in Table IV, where 
we have rounded off the 1930 km and 995 km there to 2000 km 
and 1000 km, respec t i ve 1 y . If the true wavelength is less 
than these values, the satellite track must be at an angle 
to the wave normal. The ratio of the true to the observed 
wavelength is then the cosine of the angle between the track 
and the wave normal. Horizontal wavelengths from 500 to 
2000 km were considered adequate for the purposes, and the 
plot only extends over this range. 

Figures 41 and 42 contain contours of frequency vs hor- 
izontal wavelength for the phase angles in Table IV. Thus, 
Figure 41 shows four contours with different n i trogen-oxygen 
collision frequencies V , as indicated, to obtain a nitrogen 
to oxygen transfer function phase shift of 11° as required 
in Table IV for fluctuations from the 0.003 Hz filter. Fig- 
ure 42 shows a similar set of contours for a nitrogen to 
oxygen transfer function phase shift of 21°, as required in 
Table IV for the 0.006 Hz filter. The dashed line in Figure 
41 is for the ion temperature transfer function phase angle 
of -169° given in Table IV for the 0.003 Hz filter. The 
frequency for this parameter is almost independent of dip 
angle and ion-neutral collision frequency. It is labelled 
'all T . ' for this reason. If the dominant wave system in 
the output of the 0.003 Hz filter were a plane traveling 
wave driven by the neutrals, in which the observed on-track 
wavelength were about 2000 km, as required by the 0.003 
f 1 uctuat i ons, spectra and correlation functions, the ion tem- 
perature contour would intersect the n i trogen-oxygen con- 
tours at some point. Since there is no intersection, it 
would appear from the figure that a simple plane wave cannot 
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explain the data. An intersection Mould require more than a 
tenfold increase in v Mhich is unreasonable. Further 
decrease of wavelength would not help as well. Reasonable 
tolerances for the phase shift values of Table 10 do not 
help e i ther . 


Figure 42 also contains an ion temperature contour as 
in Figure 41. It also contains some electron density con- 
tours for ion-neutral collision frequencies of 
0.44sec“ ,0 .87sec~ l and 1.75sec _1 for dip angles of 20° and 
60°. The electron density response contains a complication 
not evident in the ion temperature response. As horizontal 
wavelength is increased a root of the denominator of the 
transfer function may pass close to a complex pole of the 
determinant of the system. Depending on which side of the 
pole is passed, one may obtain a leading or lagging transfer 
function phase shift. Lower wavelengths are found to yield 
leading phase shifts. Since the desired phase shift for the 
0.006 Hz filter is 21° leadinq, one finds that the 0.87sec“* 
collision frequency and 60° dip angle contour must end 
between 1500 and 2000 km, since the phase angles are lagging 

■^ op 2000 km. The contour for 6^=60° an< ^ v j „=0 . 44sec~^ 
is lagging for wavelengths somewhat greater than 500 km, and 
is not shown for that reason. All contours for 8^=20° pro- 
vide leading phase angles. 

Note again that the ion temperature and the electron 
density contours in Figure 42 do not intersect the 
nitrogen-oxygen contours. More than a reasonable increase 
in n i trogen-oxygen collision frequency would be required 
than is acceptable to achieve this intersection. Thus, it 
would appear that a simple plane wave model does not explain 
the on-track wave for the 0.006 filter; that is, for a wave 
of about 1000 km based on fluctuations, spectra and correla- 
tion functions. Decreasing horizontal wavelength below 500 
km also does not help, as can be seen. 


In all the cases shown in Figures 41 and 42 and for 
similar computations at other temperatures, one finds that, 
if the phase shifts of Table IV are satisfied, the ampli- 
tudes of Table IV are not satisfied, the calculated values 
being too large. It is not inconceivable, though, that one 
might make a comparison for amplitudes, discarding the phase 
angles as too inaccurate. As can be seen from Figure 38, a 
given amplitude may be achieved for both a frequency that is 
higher and a frequency that is lower than that for the fre- 
quency based on the phase shift. On being consistent with 
amplitude for all parameters, one comes to the same conclu- 
sion as for phase shifts; that is, a plane travel ing wave 
cannot provide the values in Table IV. In fact, the ion 
temperature amplitude is almost always very close to unity, 
whereas the values in Table IV are larger. Graphs for 
amplitudes similar to Figures 41 and 42 will not be given 
here . 



29 


The plasma transfer functions in Figures 41 and 42 
assume that the true wave is traveling in the north-south 
direction, south in the northern hrmisphere and north in the 
southern hemisphere. It is possible that the wave normal 
has an east-west component. In this case K j n equations 
<4) and <5) is the north-sou'ch component of the true 

wavenumber vector. K j n equation (3), however, is replaced 
by the true normalized horizontal wave number. In this case 
the wavelength along the north-south direction is greater 
than the true wavelength. Computations were made for this 
case, as well, for horizontal wave number vectors from 0° to 
as much as 76 ° off the north-south direction. The results 
of these calculations are not significantly different from 
those obtained for waves traveling in the north-south direc- 
tion. No traveling wave solution is found as above. The 
results will not be given here. 
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Figures 11 to 16,18,23,26 and 27 illustrate the oscil- 
latory content of the -fluctuations o-f the various parameters 
-for data passed through various high pass -filters with cut- 
off frequencies from 0.0013 Hz to 0.013 Hz. In each of 
these examples the period of the oscillations is nearly con- 
stant. There appear to be many such periods, or 
wavelengths, as indicated by the spectra of Figure 25 and 
the outstanding waves listed in Table Ilia. 

Yet, one would not expect such near constancy in the 
period of the fluctuations over the large distance travelled 
by the satellite during the 4,125 seconds of time of Figures 
11 to 16. For one, motion is over a sphere and in both hem- 
ispheres. Also, the sources for these waves may be in both 
auroral zones and not necessarily coherent. As the satel- 
lite path enters a localized source region and moves away, 
the angle of the path relative to the wave normal of the 
waves emanating from the sources changes. For a source of 
fixed horizontal wavelength, the change in this angle with 
satellite position would change the apparent periodicity, 
or, equivalently, the apparent horizontal wavelength, which 
is not observed. The apparent constancy of the oscillations 
may imply an extended rather than a localized source region, 
which appears to be in the northern auroral region, based on 
the discussion of fluctuations in Section 3.3. However, 
sources appear present in both auroral regions with the 
southern hotter. Both could be generating waves that inter- 
fere to produce the observed pattern. These sources could 
have similar spatial and temporal structures, though not 
completely coherent. Diffuse aurora is known to have some 
symmetry in both hemispheres and is of the appropriate scale 
size to excite large scale waves <Lui et al , 1982; Akasofu, 
1974) . 


The fluctuations represent spatial waves along the path 
of the satellite with wavelengths from about 4000 km to 
about 400 km. Figures 11 to 16 demonstrate the correlations 
between the fluctuations of all four parameters that were 
studied for two different dominant wavelengths. These 
wavelengths were on the order of 2000 km in Figures 11 to 13 
and 1000 km in Figures 14 to 16. Similar fluctuations and 
correlations may be shown for other wavelengths by using 
different filters. The extent of the correlation is well 
brought out in the plots of the cross correlation functions 
in Figures 32 to 36 for the data passed through the same two 
filters as for Figures 11 through 16. Correlation of nitro- 
gen and oxygen is strong as shown in Figures 32 and 33, and 
the fluctuations appear nearly in phase, confirming the 
relationship shown in Figures 11 and 14. Ant i -corre 1 at i on 
of ion temperature and oxygen density, apparent in Figures 
12 and IS, is clear in Figures 34 and 35; however, the 
degree of correlation is not as strong as that of oxygen and 



nitrogen densities. It is interesting that the ion tempera- 
ture is somewhat better correlated with the nitrogen density 
than the oxygen density. This correlation .could possibly 
relate in some way to the chemistry and thermodynaami cs of 
the region, if not unique to this orbit. It should be 
investigated -for other orbits and their data. 

Figures 13 -for the 0.003 Hz -filter and Figure 16 -for 
the 0.006 Hz -filter, to a lesser extent, suggest that the 
electron density behaves differently in the two hemispheres, 
since the phase of the -fluctuations changes on passing from 
one region to the other. This difference may arise from the 
orientation of the magnetic field. For example, at low tem- 
poral frequencies and for horizontal wavelengths in the 
range of 500 km to 2000 km, the phase shift of the electron 
density transfer function is small and leading if the wave 
is traveling toward the equator, but large and lagging if 
traveling away from the equator. These charac ter i st i cs 
appear to be like the phase shifts observed in Figures 13 
and 16. 

Nevertheless, the data do not appear to be in accor- 
dance with a pure traveling wave, as described in Section 5. 
However, there is no doubt that the constituents, both neu- 
tral and ionized, including temperature, are coupled with 
spatial characteristics over a wide range, of scale sizes 
down to hundreds of kilometers over global regions. The use 
of 15 secono data samples precludes study of coupling for 
still smaller scale sizes. Gross et al <1980) have looked 
at these smaller scales using 1.5 second AE-C data samples 
and also found correlation between parameters. 

Since the plane wave model of Section 5. did not 
account for the char ac ter i st i cs of the data, one may expect 
that the model is too simple. Perhaps, viscosity and other 
factors that were neglected are important. Beating effects, 
over and above noise, make the analysis difficult, as well. 
Also, the transfer functions of Table IV may be of question- 
able accuracy, directly affecting the modeling analysis. A 
mixed wave system of stationary and traveling waves is pos- 
sible, which could arise from both a single or two sources, 
one for each hemisphere. In the latter case the sources 
need not be truly coherent, though of similar spatial and 
temporal charac ter i st i cs . Traveling components most cer- 
tainly exist, as is known from observations of large scale 
TID's by ground instrumentation. But ground instrumentation 
may not recognize a stationary structure. More than one 
facility may be required, distributed in some way in two 
dimensions on the ground. Also, analysis of data of many 
orbits appear to indicate that these waves with large scale 
lengths are quite common in the global structure, yet such 
waves are not commonly observed on the ground (Georges, 
1968). Experiments are needed to study these matters 
f ur ther . 
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It is difficult to analyze a system of stationary and 
traveling waves with just satellite measurements, no less in 
situ data. Ground measurements are needed to augment the 
satellite measurements by providing temporal and possibly 
directional information. The satellite provides the spatial 
distribution. The ground measurements should be made at the 
time of a satellite pass over the ground station. Because 
of the difficulty of finding coordi nated measurements in the 
past, one must arrange for such experiments in advance, if a 
suitable satellite is in orbit, and for enough time to 
insure decent and useful measurements. (Such an opportunity 
exists with the Dynamic Explorer 2 satellite.) Further 
modeling efforts are desirable, particularly for propagation 
over a sphere with a superposed dipolar magnetic field, such 
as has been started by Mayr et al <1982). 


An interesting question concerns the large number of 
waves that appear to make up the distribution, as indicated 
by the number of wavelengths in Table Ilia. Nothing like 
this has been observed with TID's measured from the ground. 
Such large numbers of waves are seen for other orbits. Of 
course, a number of waves may be present for part of the 
data span. Also, some of these waves may be viewed as the 
Fourier decomposition of a structure, or wave packet, in the 
distribution sampled by the satellite and not as individual 
waves; that is, this number of waves may be merely an indi- 
cation of the spatial complexity of the distribution result- 
ing from complex source regions. Nevertheless, each wave 
has its own identity in a linear system, as long as linear- 
ity applies. There is no guarantee that any or all are pro- 
pagating, stationary or a mixture. 

An important find is the power law decrease of spectra, 
with each parameter having its own index, though that of 
nitrogen and oxygen are close. Since the wavelength range 
is only 16.5:1 here, the fall off may not be 1 ogar i thm i cal 1 y 
linear over a larger range of wavelengths. This type of 
behavior was also observed for shorter wavelengths using 
AE-C 1.5 second data samples (Gross et al , 19S0) for an 
additional 10:1 wavelength range. Thus, it appears that 
this type of spectrum occurs for wavelengths ranging from 
thousands of kilometers to tens of kilometers, but not 
necessarily with a constant spectral index over the entire 
range. Differences in the fall off rate may arise from the 
transfer functions between major and minor species 
responses. The fall off rate itself may arise from the 
nature of the source and propagation in the intervening 
medium between source and observation point. The relation- 
ship to turbulence is of interest, and these matters deserve 
further research. Perhaps, the Van Zandt (1982) thesis for 
a universal spectrum for buoyancy waves is very revel ant. 

This study pro< ides a measure of the rms fluction level 
for the four parameters which depend on the wavelength of 
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the di sturbancee . Values are given in Table' I and II. If 
temporal -frequencies are known (even -for standing waves), 
the flux of energy and momentum may be deduced. Such deter- 
minations can only obtained from more comprehensive, coordi- 
nated measurements of satellite and ground facilities. The 
methods used here permit a study of the magnitude of fluc- 
tuations and their spectral properties under both quiet and 
disturbed conditions, as well. 

Finally, this study of AE-C data from orbit 5133 on 
January 20, 1975 is not a unique set of data. Less complete 
studies were made of data from many other orbits of that 
satellite and of data from the AE-E satellite. The results 
shown are typical. Similar fluctuations are observed, spec- 
tra are similar with peaks at particular on-track 
wavelengths and with a decreasing power law fall off. 
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Fig. 2 Electron and ion temperatures vs. time in seconds UT. 

The dotted curve is the electron temperature, the other curve 
the ion temperature. Atmosphere Explorer C UA data for 
orbit 5133, January 20, 19V5. 
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Same as Fig. 5 but for ioi. temperature. Average is 932.47°K. 
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order= 40 of oxygen and nitrogen vs. time in seconds UT. 
The solid curve is for oxygen, whereas the dotted curve is 
for nitrogen. Both of the fluctuations have been normalized 
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The aame as Fig. 11, but the dotted curve is for the ion 
temperature fluctuations. 
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Fig. 11, except the cutoff of the 
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Fig. 14, but the dotted curve is for the ion 
fluctuations. 


Electron density 



The same. *s Fig. 14, but the dotted curve is for the electron 
density fluctuations. 



Final Prediction Error (fpc) vs. MEM order for nitrogen 
density fluctuations from a high- pass filter of order = 40 and 
cutoff frequency of 0.0015 Hz. 
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Fluctuations of nitrogen density vs. time at the output of a 
high-pass filter of order = 40 and cutoff frequency of 0.0015 Hz 
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nitrogen density fluctuations from a high-pass filter of order 
= 40 and cutoff frequency of 0.0015 Hz. The order of the 
MEM is 20. The scales and numbers above the peaks of the 
spectrum are as explained for Figure 19. 
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Fig. 22 Same as Fig. 20 but for a high-pass filter cutoff of 0.003 Hz. 



luctuations of nitrogen density as Li Fig. 18 , but for a high- 
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Fig. 27 The same as Fig. 18, but for a high- pass filter cutoff 




Autoregressive model of order equal to 20 of the oxygen 
density fluctuations as expected from a high-pass filter with 
a cutoff frequency of 0.003 Hz. 
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Fig. 29 The same as Fig. 28, but for the nitrogen density fluctuations 






0.003Hz 


0 WQ,\7- ' 

0F p °°* qual 




* 







0.006Hz 


ORfCJrWL I 
OF POO,"; Qt; ll -.;vY' 



•A 











The same as Fig. 32, but for oxygen density and ion 
temperature fluctuations. 
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The same as Fig. 39 but for the ion temperature transfer 
function. 3*s and 6's on the phase shift curve are for the 
phase angles required in Table IV. See text. 
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Plot of theoretical frequency in mHz vs. horizontal wave- 
length in km to satisfy the phase shift requirements of Table 
IV for the 0.003 Hz high- pass filter. The temperature is 
830°K. The solid lines are for the nitrogen transfer function 
phase shift for the different values of the nitrogen-oxygen 
collision frequency v n . All points for the ior? temperature 
transfer function phase shift are on the dashed line. There 
are no electron density curves for this case. 
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The same as Fig. 41, except that the filter requirements 
are those for the 0.006 Hz cutoff frequency. Electron density 
contours for the transfer function phase stfift are shown for 
dip angles of 20° and 60° for propagation towards the equator. 
Ion neutral collision frequency, v^ n , values associated with 
the curves are shown. These electron density contours are 
labeled, 




